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ABSTRACT 


Techniques have been established and a computer program has been written (in the 
experimental language, FORMAC) to symbolically reduce arbitrary block diagrams for 
desired transfer functions. Symbolic solutions are determined in several forms including 
an expanded form in terms of the driving frequency and system constants. Programs are 
written to numerically evaluate the symbolic solutions for real and imaginary parts and 
magnitude ratio and phase angle. The programs have been applied to several research 
problems which include both lumped and distributed parameter systems. The latter 
forms are built into the program and are handled automatically. 
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COMPUTER PROGRAM FOR SYMBOLIC REDUCTION OF BLOCK 
DIAGRAMS USING FORMAC 
by Carl F. Lorenzo and Paul Swigert 
Lewis Research Center 

SUMMARY 

Techniques have been established and a computer program has been written (in the 
experimental language, FORMAC) to symbolically reduce arbitrary block diagrams for 
desired transfer functions. Symbolic solutions are determined in several forms including 
an expanded form in terms of the driving frequency and system constants. 

Programs are written to numerically evaluate the symbolic solutions for real and 
imaginary parts and magnitude ratio and phase angle. 

The programs have been applied to several research problems which include both 
lumped and distributed parameter systems. The distributed-parameter forms are built 
into the program and are handled automatically. 


INTRODUCTION 

The use of block diagrams has become widespread for the analysis of dynamic prob- 
lems. The classes of dynamic problems handled with this technique include controls and 
servomechanisms, various physical dynamic problems from acoustics to structural dy- 
namics, and such engineering tasks as transient heat-transfer and circuit analyses. 

In cases of block diagrams with few blocks or paths’, symbolic solution is rather 
routine. However, as the path and block complexity increases, manual symbolic solu- 
tion becomes more difficult and, indeed, prohibitive. The manual reduction of a block 
diagram is an algebraic task which is at best routine and can often be difficult and/or 
tedious. Also, the process is susceptible to human error and is time consuming. 

In the past, these more complex problems have been attacked by use of computer 
programs to allow at least numerical solutions. These programs have used matrix 
methods and require a complete matrix manipulation for each frequency (refs. 1 to 3). 



The purpose of this study is to evolve a digital computer program that can reduce any 
block diagram symbolically to obtain a transfer function of interest, that is, to generate 
expressions for the transfer function either in terms of the G's (the transfer functions 
for the blocks), or S (the Laplace operator), or a> (the driving frequency) . 

A further goal is to supply programs by which the symbolic solutions can be effi- 
ciently evaluated numerically for magnitude ratio, phase angle, and real and imaginary 
parts as functions of frequency. 

^ Such programs would allow the engineer to focus his attention on detailed evaluation 
of the model (contents of each block) and interpretation of the end results of the analysis . 

The value of the symbolic expression is that (1) it may be used in research applica- 
tions where such expressions are useful for some further analytical purpose, (2) the 
mathematical expression is a most compact manner for communicating the desired infor- 
mation (transfer function) for reports, etc. , and (3) the analytical expression allows effi- 
cient numerical evaluation for cases of particular interest. 

The present work differs from the previous efforts (ref. 1) in that symbolic solutions 
will be generated (in minimum form) using FORMAC. These symbolic solutions will be 
evaluated directly for numerical results, therein bypassing the matrix manipulation. 

This report is organized in the following manner. The first section will be general 
background information with some theoretical preliminary material to show how block 
diagrams are reduced. The next two sections will give detailed information as to how the 
general ideas are implemented into digital programs and how those programs function. 

An applications section demonstrates program results for both lumped parameter and dis- 
tributed parameter systems. A users manual is presented which gives a step-by-step use 
of the program with an example . 


GENERAL SCHEME 

The basic elements composing block diagrams are 

(1) Summers - add or subtract n-signals 

(2) Blocks - accept one signal and modify it by an S or frequency-sensitive operator 

(3) Nodes - split a signal into two or more parts 

Any linear block diagram can be composed of these elements in various combinations. 

The information contained in the block diagram when the blocks have some gener- 
alized form (G-form) is basically topological information (ref. 4). When the expressions 
for the G's are put into the diagram, the information becomes system information. 

Consider the problem of finding an arbitrary transfer function X /X, for some 
given block diagram. (Symbols are defined in appendix A.) For a sufficiently general 
block diagram, there will be several solutions: (1) the forward transfer function, (2) the 
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Figure 1. -Two-element block diagram. 


external-loop transfer function, and (3) the transfer function with respect to a given in- 
put. The forward and external-loop transfer functions are illustrated in the simple block 
diagram of figure 1. It is characterized by the equations: 


X 1 = G 1 X 2 


and 


X 2 = G 2 X 1 

To form the transfer function X^/Xg, clearly, two results are possible: 



and 



using the first and second equations, respectively. This is also true of larger, more 
complex diagrams. Generally, it is the forward transfer function (X^/^ = G^) that is of 
interest in the solution of controls problems. When a block diagram is reduced manually, 
the problem does not usually occur because the engineer is discriminatory in his selec- 
tion of equations (blocks or summers) used to reduce the system; that is, he takes into 
account the direction of signal flow. It is undesirable to account for signal flow direction 
in the computer program. Therefore, it is best to avoid this problem by forming only 
transfer functions involving an input, as opposed to those containing two internal signals 
in the block diagram . 

The effect of an input on the solution is related to the question of forward and 
external-loop transfer functions. To illustrate, the block diagram of figure 2 has two 
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Figure 2. - Two-input block diagram. 


inputs. As a result of this, there will be two transfer function relations for X 2 /X 4 , for 
example. Consider the system equations: 


( 1 ) 

( 2 ) 

(3) 

(4) 


x 5 i + Gi g 2 

also, 

X 4 _ ~ G 2 

X 5 1 + G 1 G 2 

Hence, 


Xj - Gl x 4 = 0 


X 3 - G 2 X 2 = 0 


X 6 - X 3 " X 4 = 0 


X 5 + X l - X 2 = 0 


It is easily shown that for Xg as input (Xg = 0) 


It 


1 


x 2 1 + 0^2 



1 + GjG 2 

can be further shown that for Xg 



for X 5 


as input (Xg = 0) 


as an input 
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and 


Hence, 


G, 


X 6 1 + G 1 G 2 


X, 


X, 


1 + G i G 2 


1 + G 1 °2 

= — — = for Xg as an input 

1 + G 1 G 2 

It is concluded, therefore, that the transfer function between two arbitrary points 
(neither an input) in a block diagram requires the specification of an input as reference; 
that is, total specification of a transfer function should be X a /X b with X c as an input 
(which implies that all other inputs are zero) . Or how does X & respond to X^ when 
X c (only) is stimulated? It is interesting to note that, when both inputs Xg = Xg = 0, 
equations (1) and (4) yield X 2 /X 4 = G 4 and equations (2) and (3) yield X 2 /X 4 = -l/Gg. 

These, by analogy to the previous example, would be considered the forward and 
external-loop transfer functions. This also indicates that, if all the equations are not 
required for a reduction, the solution is probably not unique. Were it not for this fact, 
it would be difficult to determine whether the solutions for large complex diagrams were 
unique . 

The problem of nonunique solutions is easily remedied in the programs which follow 
by only forming transfer functions with respect to an input. Experience indicates that this 
results in a unique solution. Ratios of the proper results form the arbitrary transfer 
functions; that is, 



X„ 
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m 


X„ 


X v 


X. 
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for input X m 
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Numerical manipulations 



Figure 3. - Simplified flow chart for block-diagram-reduction technique. 


It is further noted that, if minimum form solutions are attained for X a /X m and X b /X m , 
the denominators (characteristic parts) will generally be identical. Hence, only the ratio 
of the numerators are needed for the desired result. This can be validated by applica- 
tions of Cramer's rule. 

A simplified flow chart for the process of reduction of a block diagram (to a desired 
transfer function) is presented in figure 3. The flow chart is composed of two basic 
parts: (1) those processes involving algebraic manipulations and (2) those processes in- 
volving numerical manipulations. 


Algebraic Manipulations 

The algebraic manipulations, of course, precede the numerical. The first block 
(fig. 3) indicates the information required for the algebraic (or symbolic) reduction. 
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r 


r These are (1) the equations corresponding to each element of the block diagram (one 
:^r- equation per element), (2) the transfer function sought from the reduction, and (3) the 
: contents of the blocks of the diagram, that is, the equations expressing the G’s in 
terms of the Laplace variable S and the various time constants and natural frequencies, 
etc. These data are sufficient to allow reduction to a symbolic form. 

Before the actual algebriac reduction can take place, the order in which the variables 
are to be eliminated and the equations to be used for such eliminations must be deter- 
mined. For a manual reduction this is usually done by observation or intuition. The 
technique of ordering will be discussed more fully in appendix D. 

Having a desirable order of elimination, the reduction takes place in the following 
manner. The first variable to be eliminated is solved for in the equation indicated by the 
order determination. This equation is then substituted into all other equations of the set 
containing that variable, thereby eliminating it. It is important that the variable be sub- 
stituted in all equations; otherwise, the possibility of looping occurs. (Looping is re- 
peated substitution of variables without achieving elimination.) This substitution process 
is continued for all the remaining variables except the two variables involved in the de- 
sired transfer function. This leaves a single equation in the two variables which is now 
manipulated to form a ratio only involving the G's. For example: 

— = ± (5) 

X b 1 + G 1 G 2 

This is the so-called G-form of the transfer functions solution. The G-functions can 
always be expressed as a ratio, that is, G = N /D . This ratio is now substituted into 

a d, a 

equation (5) which results in 





( 6 ) 


Solving for a simple ratio for Xg/X^ gives 


D 1 D 2 


X b D 1 D 2 + N 1 N 2 


( 7 ) 


for the example taken. In this form, the system information is easily introduced into the 
result. Assume that 
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and 


TjS + 1 
t 2 S + 1 




>- 



t 3 S + 1 
t 4 S 2 + 1 


J 


( 8 ) 


The fourth block (flow chart, fig. 3) indicates that equation (8) is to be substituted in 
equation (7), thus, 


X a _ + 1 )( T 4 S2 + *) ,, 

X » ( t 2 S + ‘)( T 4 S2 + ») + (V + 1 )( t 3 S + >) 

This is the S-form solution for the transfer function (unexpanded). From this form, 
either numerical evaluation is possible or further symbolic reduction can be done. For 
the symbolic result, S is now replaced by io> (a step which is valid for wide set of con- 
ditions) . (See ref. 5.) Hence, 


X„ 


(r 2 iw + l)(r 4 i 


(t„ io> + 3 ) (x^i 2 ^ 2 + l) 


X b (t g io> + l)(r 4 i 2 a , 2 + l) + ( T ^iw + l^Tgioj + l) 


2 

Now, i = -1, thus 


X„ 


(r 2 ia» + l)(- 


2 i 

T„W + 1 


' r A ) 


^b + 1 ^-t 4 oj 2 + + ^io; + l^Tgicc + 1^ 


and 



(l - r 4 o; 2 ) + i(- 


r 2 r 4 cv 


r 2 cv) 


2 ~ ( t 1 t 3 + t 4) oj2 +i ( T l +T 2 + ^ " t 2 t 4 w3 
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This has been called the complex rational form in the sixth block of the flow chart 
(fig. 3), that is, the form 


= A ± iB (10) 

X b C + iD 

where A, B, C, and D are functions of the excitation frequency and system parameters, 
t ’ s , o> n 's, K's, etc. This form is very important in many research studies. From this 
form, the magnitude ratio and phase angle follow directly because 


Re 




A C + BD 
C 2 + D 2 


CB - AD 
C 2 + D 2 


(U) 

( 12 ) 


therefore, 


Phase angle = tan 


-1 


/ CB - AD \ 
\AC + BD/ 


(13) 


also, 


Magnitude of 





/a 2 I B 2 

y c 2 + D 2 


(14) 


Numerical Manipulations 

For the numerical manipulations, when the path starting with the S-form is con- 
sidered, it is merely necessary to substitute the numerical values for the frequency and 
the system parameters then to use complex arithmetic to obtain the numerical equivalent 
of equation (10). From this, the quantities of equations (11) to (14) are readily evaluated. 

The path from the complex rational form is similar except that the substitutions are 
made into equation (10) directly. 

A further note on the importance of the order used in the reductions: if the proper 
order is not used in the reduction scheme presented, it is possible to introduce a common 
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factor into both the numerator and denominator (of eq. (9), e. g.). Now, although this re- 
sult is still technically correct, this could potentially introduce hundreds or thousands of 
extra terms in the symbolic solution. This is important in the computer program which 
follows because it seriously affects storage capabilities and, hence, problem size capa- 
bility . 


Problem Types 


The type of problems that can be handled using these techniques is limited to linear 
systems. Hence, the contents of the blocks (the G's) can only be functions of S and 
constants; that is, time varying constants (system characteristics) are not permissible. 

For the so-called lumped parameter systems, the G’s are sums and products of 
KS terms and constant terms. 

For distributed parameter systems, namely, systems where the blocks are solutions 
to the wave, diffusion, or beam equation or where they are dead times, the G’s will 
generally occur in the following forms (one dimensional) : 


sinh (d}fi 


sinh d ? as + bs + c 


(-v: 


cosh (d T as + bs + c 
sin (d Vas 2 + bs + c ) 




COS ( d T as + bs + c 

V ' 

d Vas^ + bs + c 


exp 


Vc - 2 


as + bs + c 


(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 


Here again, this is true provided the partial differential equation involved is not of the 
time varying coefficient type. 

The preceding forms are built into the programs described in the next section. The 
following identities are useful when these forms are involved in either a manual reduction 
or in the programs. The program symbolic output is RSINH (real part of sinh), ISINH 
(imaginary part of sinh), etc. , when z = x + iy. 
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sinh z = sinh x cos y + i cosh x sin y (21) 

cosh z = cosh x cos y + i sinh x sin y (22) 

sin z = sin x cosh y + i cos x sinh y (23) 

cos z = cos x cosh y - i sin x sinh y (24) 

e z = e x cos y + i e x sin y (25) 


These equations can be found in reference 6. In the event of a square-root argument, 
the following equations are helpful: 

fz = Vx + iy = u + iv (26) 


If x>0 


If x<0 



(27) 


u = y/2v ^ 


v = (sgn y) 



+ -|x + iy| 
2 


> 


J 


(28) 


The forms (21) to (25) result when the various partial differential equations are solved in 
rectangular coordinate systems. When this is not the case, other forms can be intro- 
duced, for example, Bessel functions and Legendre functions, etc. When this is the 
case, it becomes desirable to have some arbitrary unassigned functions which can repre- 
sent such a solution. It is only necessary that expressions (similar to (21) to (25)) can be 
written for the real and imaginary parts. It is interesting to note that the sinh, cosh, 
sin, and cos will not have square-root arguments unless damped versions of the original 
partial differential equations (wave or beam) are solved. Whereas, the square-root form 


is usually required for the diffusion equation solutions, that is, sinh 
exp (dVi). (See ref. 7 for typical results.) 



and 
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IMPLEMENTATION 


The complete solution of the block-diagram reduction is obtained through the use of 
two separate computer programs: the algebriac solution and the numerical solution. 

The algebriac computer program required the use of the FORmula MAnipulation Compiler 
(FORMAC) (ref. 8). The FORMAC language is an experimental extension of FORTRAN 
IV for the IBM 7094 computer. This extension provides the capability for doing algebriac 
manipulation by using symbol manipulation. The numerical solution of the FORMAC ex- 
pressions is performed by a program written in FORTRAN IV. Appendix B presents a 
User’s Manual for these programs and appendix C gives listings, flow charts, and de- 
scriptions of the programs . 

Two separate programs were necessary for two reasons: (1) FORMAC is an inter- 
pretive language that requires relatively large computer execution times for repetitive 
numerical evaluations of algebriac expressions and (2) computer storage must be saved 
because the FORMAC program requires large amounts of storage to form the algebriac 
expressions. These two considerations made the separate programs necessary even 
though the user was given the inconvenience of making two computer runs for the com- 
plete solution. 

The solution obtained from the FORMAC program is communicated to the FORTRAN 
program by equations punched on cards. These equations are punched by the FORMAC 
program in a form that can be compiled by the FORTRAN IV compiler. The cards are 
simply inserted in the proper place in a FORTRAN subroutine. 


FORMAC Program (REDUCE) 

Basically, the FORMAC program consists of four major steps: (1) reduction of the 
system of equations (of the block diagram) to one equation containing the two variables of 
the desired transfer function, (2) solution of this equation for the transfer function in 
terms of the G-functions, (3) substitution of the system information for the G-functions, 
and (4) expansion of this equation and collection of terms of the real and imaginary parts. 

The reduction proceeds by choosing a variable to eliminate, solving for that variable 
in one of the equations, and substituting the solution for that variable in the remaining 
equations. Prior to the actual reduction, the inputs to the block diagram, supplied by the 
user, are replaced by zero. If an input occurs in the desired transfer function, the input 
is not made zero. With this method, the number of equations and variables is reduced by 
one at each step. The variables that appear in the transfer function are, of course, 
never eliminated. The final equation contains only the variables of the transfer function 
and the G-functions. 
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TABLE I. - STEP USED BY FORMAC PROGRAM TO SOLVE HEAD RESPONSE 


BLOCK DIAGRAM FOR TRANSFER FUNCTION Xj/Xg 



This reduction method is depicted, using the example of appendix B, in table I. 

Here, the transfer function desired is Xj/Xg. The equations to be solved are listed at 
step 0. As the reduction proceeds through each step, the number of variables and equa- 
tions is reduced by one until only one equation remains. One can see that equation (2), 
as it appears in step 4, is indeed the desired solution. 

This method will always give a correct result. However, the answer may be of a 
form that is overly complicated; that is, the numerator and denominator may contain like 
factors. Because these like factors are difficult to detect in FORMAC, it was necessary 
to choose the variables to eliminate and the equations used to solve for these variables 
in some specific order. The order desired is the one that yielded transfer functions 
which had no common factor in the numerator and denominator. In other words, the 
minimum form of the solution is desired. The method for determining this order is de- 
scribed in appendix D. 

The FORMAC program has certain limitations because of programming considera- 
tions and computer memory size. There is an upper limit to the number of equations, 
the number of variables, and the number of G-functions that it can handle. These limits 
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are 30, 63, and 20, respectively. There are also limitations on the number and type of 
functions and variables that may be used to describe the system information. All limi- 
tations are described in the user’s manual (appendix B). The limitations are due partly 
to the method employed to order the equations for solution and partly to conserve memory 
locations. It is felt that the limits are large enough so that the user is not unduly re- 
stricted in the size of problem that may be solved. However, many of the limits may be 
increased by changing dimension statements in the program. Care must be taken in 
changing the dimensions so that too much storage space is not removed from the space 
needed to store the FORMAC expressions. 

Because FORMAC stores expressions in the computer memory locations not used by 
the programs, it is desirable to keep the unused storage as large as possible. Subrou- 
tines referenced by FORMAC programs but not executed in this program are not loaded 
at execution time. This is accomplished by the use of a dummy subroutine with entry 
point names the same as those subroutines that would be loaded but not executed. 

In order to expand the transfer function into real and imaginary parts, it was neces- 
sary to represent the allowed system functions symbolically as complex numbers. This 
is accomplished in the program by prefixing the function names with R and I, where 
the names prefixed with R stand for the real part of the complex number and I, the 
imaginary part. The term S is considered as a pure imaginary number ico. 

The expansion of the transfer function into real and imaginary parts may not always 
be accomplished by the program. The limiting factor here is memory size. Because of 
this fact, provision is made in the FORTRAN program to accept transfer functions in 
either the S form or the complex rational form. 

The FORMAC program uses the following function subprograms that are not supplied 
with either FORMAC or FORTRAN. These subprograms are not given in this report and 
must be supplied by the user. 


Function 


Description 


ALS(N, X) 
IARS(N, X) 
LGR(N, X) 
AND(X1,X2) 
OR(Xl,X2) 


Accumulator left shift of X, N binary places 
Accumulator right shift of X, N binary places 
Logical left shift of X, N binary places 
Logical intersection of XI and X2 
Logical union of XI and X2 
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FORMAC Input/Output 


The input to the FORMAC program consists of two parts . The first is the transfer 
function desired and the system of equations from the block diagram. The second part is 
the system information expressions. After the equations are read, they are reduced to 
the G-form and output printed and punched. The G-form is then transformed into the 
N/D-form and output. A logical variable is read to determine whether the user wants 
system information substituted for the N/D’s. If substitution is desired, this informa- 
tion is read and the substitution and the output of the S-form takes place. Another logi- 
cal variable is read and tested before expansion to the complex rational form is per- 
formed. A more detailed description of the input and output is given in the users manual 
(appendix B) . 

For the equations punched by the FORMAC program to be acceptable to the FORTRAN 
IV compiler, the floating point powers and the trailing $ had to be removed. This is 
done by the program just prior to the punching of equations. 

The execution time for the FORMAC program is difficult to predict. However, exe- 
cution times cited in the next section are typical. All execution times mentioned in this 
report are for the IBM 7094n~7044 direct-couple system. 

FORTRAN Program (EVAL) 

The FORTRAN program can accept punched output from the FORMAC program in 
one of two forms: either the unexpanded (or S-form), or the expanded (or complex ra- 
tional form). The FORTRAN program is designed to handle all the different types of sys- 
tem information acceptable to the FORMAC program except the arbitrary function. 

To be able to handle the arbitrary function, the user must supply a complex function 
subprogram F(I) that evaluates the desired function. The subprogram that defines the 
other functions (appendix C) may be used as an example of how this is to be done. 


FORTRAN Input/Output 

The equations punched by the FORMAC program are inserted into one of two- subrou- 
tines, depending on the form of the equation. This subroutine must be compiled and exe- 
cuted with the main program and the function evaluation subprogram. The numerical in- 
put is read by the main program, and the subroutine containing the FORMAC expressions 
is executed. The numerical input is determined by the type of functions used in the sys- 
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tem information expressions. A more detailed discussion of this may be found in appen- 
dix B. 

Output from the FORTRAN program consists of all numerical input and the follow- 
ing data from the transfer function at each frequency specified by the user: 

(1) Values of the real and imaginary parts of the numerator and denominator nor- 

malized such that either the real or imaginary part of the numerator is one. 

The normalization is necessary because some terms may involve large values 
of frequency raised to powers. 

(2) The real and imaginary parts 

(3) The complex absolute value (magnitude) 

(4) The complex argument (phase angle) 

The execution time of the FORTRAN program depends on the complexity of the trans- 
fer function and is therefore difficult to predict. The execution time, however, is some- 
what larger for the S-form than the complex rational form. 


APPLICATIONS 

Several applications have been chosen to demonstrate various aspects and poten- 
tialities of the block-diagram reducing programs. The first example is that of an 
interacting-control system. It is a lumped parameter example in which the forms of the 
solutions are of some interest. A second example, representing a lumped mass, spring- 
damper system, is presented. Here, the symbolic solution is of particular interest. 
Finally, a distributed parameter is presented in which the numerical results are of pri- 
mary interest. 


Application 1: Interacting-Control System 

A block diagram of an interacting- control system is presented in figure 4(a). Al- 
though the values for the G’s are not specified here, the block diagram might be repre- 
sentative of the interacting-control problem for jet engines (e. g. , attempting to control 
engine speed and turbine -outlet temperature simultaneously). The algebriac information 
required, the program, and the solutions desired are indicated in figure 4(b) . Because 
the G's are not specified, the solutions cannot give system information, but it does give 
the topological information associated with the form of the diagram . 

Because there are eight components in the block diagram (summation points and 
blocks) there are eight equations required to define the system. (Nodes are not consid- 
ered components because the equation X = X„ is not needed by the program.) 

ct a 
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(a) Block diagram. 


Xi - X 3 - X 5 = 0 

Xy - GjXij = 0 

x 2 - x 4 - x 6 = 0 

Xg ~ G2X5 = 0 

x? + x 1Q - x 3 - 0 

x 10 ' G 3 x 6 = ( 

x 8 + x 9 - x 4 = 0 

x 8 ' G 4 x 6 = 0 

Block-diagram inputs: 

Xj and x 2 


x, x, x. x« 

Solutions desired: -2, - 2 , A — 

X 1 x 2 X 1 x 2 

(b) Algebraic information required. 

Figure 4. - Interacting control system. 


The solution for the transfer functions are as follows: 


X 4 


G 2 

X l" 

G 1 

+ - GgGg + G^ + 1 

X 3 


G 1 + G 1 G 4 " G 2 G 3 

X 1 

G 1 

+ G j G 4 - G 2 Gg + G^ + 1 

X 4 


G 1 G 4 " G 2 G 3 + G 4 

X 2 

G 1 

+ G j G^ - G 2 Gg + G^ + 1 

X 3 _ 


G 3 

<N 

X 

G 1 

+ G^G^ - G 2 Gg + G^ + 1 
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These are the results as generated by the computer. Should the user be interested in the 
transfer function X 4 /X 3 , it follows that, for Xj as the reference input, 


Go 


X 3 G 1 + G 1 G 4 " G 2 G 3 


and, for X 2 as the reference input, 



G 1 G 4 


G 2 G 3 + G 4 
G 3 


clearly not the same result. For this reason, the user should form the transfer func- 
tions as indicated herein so that the reference input will be known. 

The total computer execution time required for these four solutions was 0. 30 minute. 
It will be noted also that the solutions are in minimum form. Also, the equations are not 
put into the program in any particular order, but the order in which they are used is de- 
termined in the program. 


Application 2: Mass-Spring-Damper Dynamics 

This particular application was taken from a study in which the symbolic solution was 
of primary interest for use in further analysis. (See ref. 9.) The object of that research 
was to relate the individual dampers of n-mass systems (n = 3, fig. 5(a)) to a modal 
equivalent damper. To accomplish this required the symbolic solution in the complex 
rational form. 

The block diagram for the system is presented in figure 5(b) . The values for the 
G's are given in the blocks and variable names used in the program are parenthetically 
noted after each variable . 

For this problem, the following symbols apply: 

B damper constant, (sec)(lb)/ft; (sec)(N)/m 
F force, lb; N 
K spring constant, lb/ft; N/m 
M mass, (sec^)(lb)/ft; (sec^)(N)/m 
y displacement, ft; m 
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(a) Physical system. 


(b! Block diagram. 


g 1 x 8 ‘ X 1 * 0 
G 2 x 3 - *8 ” 0 

x 2 - x i - x 3 “ 0 

G 3 x 7 ~ x 2 = 0 

x 10 - x 8 ' x 7 = 0 
Block diagram input: x 6 


g 4 x 5 ‘ x 10 = 0 
x 4 - x 2 - x 5 = ° 

g 5 x 9 - x^j • 0 

X 11 + x 6 " x 10 x 9 = 0 
g 6 x 4 ' X 11 ' 0 


Solution required: 

x 6 

(c) Algebraic information required. 
Figure 5. - Three-mass spring -damper. 


Specifically required in the study were the coefficients of the various powers of at 
in the characteristic equation (denominator of the complex rational form) . 

The input equations required by the program are shown in figure 5(c). This is a 
10-component (10 equations) system and the solution required is for Xj/F = (X 1 )/(X g ). 
The solutions as generated by the program are as follows: 

The G-form of the transfer functions. 


V x 6 * G l G 2 G 3 G 4 G S /( 


G 1 G 2 + G 1 G 2 G 3 G 4 “ 


G 1 G 2 G 3 G 4 G 5 G 6 + G 1 G 2 G 4 G 5 ' 


G 1 G 2 G 5 G 6 


+ G 2 G 3 + G 2 G 3 G 4 G 5 - G 2 G 3 G 5 G 6 + G 3 G 4 " 


G 3 G 4 G 5 G 6 + G 4 G 5 “ 


G 5 G 6 + 1 ’ 


0 ) 
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After program substitution of the G equations: 


Xj/Xg = (K 1 + SZj)^ + SZ 2 ) / [-(K t + SZj)^ + SZ 2 )(-K 3 - SZg) + (Kj + SZj) 

x (K 2 + SZ 2 )S 2 T 1 + (K x + SZj)^ + SZ 2 )S 2 T 2 + (K 1 + SZ 1 )(K 2 + SZ 2 )S 2 T 3 
- (K t + SZjH-Kg - SZgjS 2 ^ - (K x + SZjJf-Kg - SZ 3 )S 2 T 2 + (Kj + SZ^S 4 
X TjTg + (Kj + SZ 1 )S 4 T 2 T 3 - (K 2 + SZ 2 )(-K 3 - SZ 3 )S 2 T 1 + (K 2 + SZ 2 )S 4 
X Tl T 2 + (K 2 + SZ 2 )S 4 Tl T 3 - (-Kg - SZ 3 )S 4 Tl T 2 + S 6 Tl T 2 T 3 ] 


After program substitution of S = ico and expansion to the complex rational form: 

x f/ x 6 = K 1 K 2 “ c ° 2z i z 2 + i ( K i Z 2 + K 2 Z 1^/ { K 1 K 2 K 3 " T 1 T 2 T 3 W + ^ K 1 T 1 T 3 + K 1 T 2 T 3 

+ K 2 T lT2 + KgTjTg + KgTjTg + + TjZ-^Zg + TjZgZg + TgZjZg 

+ T 2 Z 1 Z 3 + T 3 Z 1 Z 2>“ 4 + <- K X K 2 T l - K 1 K 2 T 2 ‘ K 1 K 2 T 3 " K 1 K 3 T 1 ' K 1 K 3 T 2 

- KjZjZj - KjKjTj - KjZjZg - K 3 ZjZ 2 )w 2 + i[( T l T 2 Z 2 + T 1 T 2 Z 3 + T 1 T 3 Z 1 


+ T 1 T 3 Z 2 + T 2 T 3 Z 1^ w + ^" K 1 T 1 Z 2 " K 1 T 1 Z 3 ' K 1 T 2 Z 2 " K 1 T 2 Z 3 K 1 T 3 Z 2 

- KgTjZj - KgTj^Zg - K 2 T 2 Z 1 - KgTgZj - KgTjZj - KgTjZg - KgTgZj 

- Z 1 Z 2 Zg)C 0 3 + (Kj^Zg + K^KgZg + KgKgZ^cJj^ 

The actual program output form corresponding to this equation is shown in table II. 

The coefficients generated here, with T n = M n and B n = Z n , compare exactly with 
those generated manually in reference 9. Thus, the tedious hours required to the various 
cases presented in that study can now be avoided by use of this program. 

The symbolic form of the solution is fairly complex; however, parametric studies 
using the solution are now easily made. In this case, for example, numerical values for 
the B’s and the K's can be substituted in the complex rational form equation. This 
substitution would leave an expression involving only the Z's (or damping terms). From 


20 


TABLE n. - PROGRAM OUTPUT FORM 


THE NU MER AT OR 

« KC1)*S*Z( 1) )*(M 2I+S*Z* 21) 

THE i ))*(K( 2MS*Z( 2»l*«— KC3»-S*2€3» )♦ (K(1)+S*Z(1 ) ) * ( K(2 I +S* 

Z(2I >*S**2 *T ( 1 ) ♦( K( 1I*S*Z( 1)I*(K( 2I+S*Z(2> !*S**2 *T (2 ) «■ ( K (1 )+S* 
Z ( 1))*( K( 2 )+S*Z( 2 ) )*S**2 *T( 3)-(K( 1>+S*Z(U ) * (-K (3 1 -S*Z (3 1 )*S** 

2 *T( 1|-(K( 1)+S*Z( 1JI*(-M3)-S*Z(3I )*S**2 *T ( 2) ♦ I Ml ) + S*Z( 1 1) *S 

**4 *T( 1 ) *T ( 3 1 ♦( K( 1)*S*Z( 1I»*S**4 * T (2 1* T { 31- ( K (2 >+ S*Z (2 )) * (-K(3 

)-S*Z( 31 l*S**2 *T( I>+(K(2!+S*Z(2)!*S**4 * T ( 1! *T (2 ) + (K (2 J *S*Z (2 )) 

*S**4 *T( ll*T(3)-(-K( 3)-S*Z(3))*S**4 * T ( U* T( 2 ) + S**6 *T(1)*T(2> 

*T(3I 


THE NUMERATOR (REAL I 


K(1I*K( 2) 


-W**2 *Z ( 1 ) *Z ( 2 ) 

THE NUMERATOR (IMAGINARY! 

( K ( i I *Z ( 2 ) +K ( 2 l*Z ( 1I)*W 
THE DENOMINATOR (REAL! 

Kill *K( 2 l*K( 3) 


— T ( 1)*T( 2I*T(3I*W**6 


(K( 1 ) *T ( 1 ) *T ( 3 I-nKI 1 )*T( 2 )*T( 3!*K(2)*T(1)*T(2I*K(2)*T(1 ) *T (3)*K(3I* 
T(1)*T( 2 l+T( 1)*Z( 1)*Z( 21 *T( 1 )*Z( 1I*Z{ 31+ T(1)*Z(2)*Z (31* T (2! *Z(1!*Z 
( 2 ) +T ( 2 I *Z ( 1 )*Z( 3 1 4T( 3 )*Z ( 1 l*Z4 2 ) )* W**4 


(-K( 1)*K( 2)*T( 1 )-K( 1 )*K( 2)*Tl 2 HK( 1I*K(2I*T(3)-K(I)*K(3» *T(1 )-K(ll 
*K(3)*T( 2 !— K( 1)*Z( 2 )*Z( 3)-K( 2 !*K ( 3 1* T( 1 l-K ( 2 ) * Z ( 1 ) * Z (3 1 -K (3 ) *1 ( 1 > * 
Z ( 2 I ) *W **2 


THE DENOMINATOR (IMAGINARY) 


(T( 1)*T( 2)*Z( 2)+T( 1)*T( 2)*Z( 3)+T( 11* T( 3)*Z(1 >4 T ( 1 1 *T (3 ) *Z (2 ) 4T (2 ) * 
T(3)*Z( II )*W**5 


( — K( 1 l*T ( 1I*Z( 2I-K( I)*T( 1J*Z(3)-K(1!*T(2!*Z(2)-KU)*T(2)*Z(3)-K(1) 
*T(3)*Z( 2)-K( 2 )*T( 1 )*Z ( 1I-K( 2 l*T ( 1 )* Z ( 3 »-K ( 2 1 *T (2 I *Z (1 ) -K (2 > *T (3 I * 
Z ( 1 1 — K( 3 I *T( 1)*Z( 1)-K( 3)*T( 1 l*Z ( 2 l-K ( 3 1 * T( 2 1 * Z ( 1 |-Z (1 1 *Z (2) *Z (3 > I * 
W**3 


(K(1)*K( 2)*Z(3)*K( 1 )*K(3 1*Z( 2 1 *K ( 2 )* K( 3)*Z( 1 ) )* H 



V 1 

I' 

this expression, the effect of the individual dampers on the response of interest could be 
assessed. 

The computer time required to execute this study was 0. 39 minute, the user prepa- 
ration (for the program) was minimal because it merely required listing the equations 
(the G-expressions) and transfer functions sought and placement of appropriate punched 
cards in the program. 

Also, with the numerical evaluation program, frequency response for cases of 
special interest (M's, K's, and B's) are quickly generated. 

Application 3: Hydraulic-Line-Dynamics Study 

This study is an example of a distributed parametric study. For this study, only the 
numerical results were of interest. For the program, of course, it is necessary to gen- 
erate the symbolic answers at least to the S-form to obtain the numerical answers. 

The following symbol list applies to this application: 

C n line capacitance per unit length, in. ; m 

2 2 

C compliance, in. ; m 

m 

p o O Q 

L n line inductance per unit length, sec /in. ; sec /m 
l line length, in . ; m 

m index number for termination impedance 

n index number for characteristic impedance 

2 2 

P pressure, lb/in. ; N/m 

3 3 

R line resistance per unit length, sec/in. ; sec/m 

2 2 

R m resistance, sec/in. ; sec/m 

W weight flow, lb/sec; N/sec 

Z Qn characteristic impedance equal to l/S V V C nP + <V C n> 2 * * S ] , sec/in. sec/m ^ 

2 2 

Z termination impedance equal to R„/(R rri C S + 1), sec/in. ; sec/m 

m in ah in 

7 n propagation constant equal to i/in- ; Vm 

The physical system is shown schematically in figure 6(a). The system is composed 
of a main feed line OJ feeding three subsidiary ducts, JB, JC, and JD. The dynamic 
problem is to determine the response of the end pressures in the subsidiary ducts in re- 
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*2 ~ ^ 1^1 = ® 
x 3 ' G 2 x 20 ° 0 
x 19 ’ G 3 X 4 = G 
' x 18 ' G 4 x 17 = 0 
x 5 " g 5 x 14 = 0 
x 6 ’ G 6 x 4 = 0 


x 7 * G 7 x 15 ‘ 0 
x 8 ‘ g 8 x 4 = 0 
x 9 " g 9 x 16 = 0 
x 10 " g 10 x 4 “ 0 
x 14 ' G ll x ll = 0 
x 15 ' G 12 x 12 = 0 


Block diagram input: Xj 
P X 

Solution required: _£=J! 

Po X 1 

Auxiliary information: 


S 2 + R n C n S 


n = y*-n G n 


' m RmV + l 


(c) Algebraic information required. 


x 16 ' g 13 x 13 1 0 


x 2 -x 3 -x 4 -0 
x 19 + x 18 " x 20 = 0 


x 6 - x 5 - X 11 = 0 


x 8 - x? - x 12 - 0 

x 10 ' x 9 ' x 13 = 0 
x 14 + x 15 + x 16 ' x 17 = 0 


Line 

Length 

Line inductance per 
unit lenath 

Line capacitance per 
unit lenath. 

Line resistance 
Der unit lenath. 

Resistance, 

R m 

Resistance 
times com- 




I 


r 



R 





pliance 


in. 

cm 

L n 


n 


"n 







sec^/in. 3 

sec 2 /m 3 

in. 

m 

sec/in. 3 

sec/m 3 

sec/in. 2 

sec/m 2 

K irrnv 

sec 


197 

5.0 

2. 73x10 ' 5 

1.66 

12. 38x10 " 5 

31. 44x10 " 7 

0 

0 




A 




B 

72 

1.83 

13. 14 

8. 018 

.94 

2.388 





8.48 

13xl44xl0 3 

0.0697 

C 

75 

1. 90 

6. 71 

4.094 

1. 78 

4.521 





2.61 

4.045 

.0212 

D 

32 

.813 

6. 71 

4.094 

1. 78 

4.521 





2.61 

4.045 

.0212 


(d) Numerical information required. 
Figure 6. - Concluded. 


sponse to pressure disturbance at the far end of the feeder duct, for example, to deter- 
mine the frequency response of P C /P Q . 

The system block diagram is shown in figure 6(b). The contents of the blocks, be- 
cause each line is a distributed system, are hyperbolic functions of square -root argu- 
ments. The form of these G's represents four -terminal network solutions for the wave 
equation for each section of line. The hyperbolic functions are multiplied by impedances 
representing the line terminations. These impedance forms together with the equations 
for the block diagram are given in figure 6. 
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Using the above information, a computer solution for the S-form transfer function 
was generated in 1.10 minutes of computing time . (The solution requires two pages of 
printed output for the symbolic answers which are not presented herein.) A deck of 
punched cards is also output which contains the symbolic solution. 

For the numerical values of the line characteristics presented in figure 6(d), the 
frequency response was numerically evaluated for a range of frequencies of 0 to 50 hertz 
in 1 -hertz steps. The amplitude ratio and phase -angle results are plotted in figure 7. 
The execution time required to generate these data was 0. 50 minute. 


Other Applications 

The applications stemming from having a transfer function readily available in sym- 
bolic form are many. For example, if one is dealing strictly in lumped parameter sys- 
tems, a program could be evolved which would evaluate system stability quickly for a 
large number of numerical cases by application of Routh's criterion to the characteristic 
equation (which are now available in symbolic form directly from the program) . 

System stability can also be determined using the slope of the amplitude ratio plot as 
it goes through the zero dB line, or by applying Nyquist’s criterion to the real and imag- 
inary parts. 

Nonlinearities can be introduced into the general analysis by use of describing func- 
tion and small perturbation techniques, both of which give linear approximations for the 
block transfer functions. (See ref. 10.) 


CONCLUDING REMARKS 

Techniques were established and a program was written to reduce arbitrary block 
diagrams to obtain symbolic expressions for any transfer function with respect to an in- 
put. Further evaluation programs are presented which allow these symbolic solutions to 
be evaluated for a desired set of frequencies and numerical values of system constants. 
These evaluations yield magnitude ratio, phase angle, and real and imaginary parts of 
the transfer function. 

The advantages of the present study over previous techniques are as follows: (1) the 
results are given in symbolic form (3 forms: G, S, and complex rational) and (2) ar- 
rangement of the input data in the form of arrays and the matrix manipulation for each 
frequency point are not required. (The symbolic solution need be generated only once 
and then reevaluated for each frequency.) 
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The most serious limitation of the work is the computer storage required. The 
FORMAC routine itself requires a large amount of storage, and, as the systems size in- 
creases, the transfer function equations grow very large and, hence, the storage capa- 
bilities are affected. Specifically, as written, the program can handle a system of 
30 equations, containing as many as 63 variables. 

The computing time required is modest even for large systems. Also, the user 
preparation time is minimal. 

The program is capable of handling block diagrams for any linear constant coefficient 
system. The system can be either a lumped parameter or a distributed parameter sys- 
tem so long as an expression for the real and imaginary parts is known. Furthermore, 
the ability to handle distributed parameter block in rectangular coordinate form is built 
into the program. 

Finally, a criterion is presented which allows ordering of the equations and variables 
to be eliminated so that a minimum form solution can be achieved. This technique has 
given minimum form solutions for many trial cases. 

The technique involved for the algebraic reduction is simple and handles the problem 
of looping. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 26, 1968, 

180-31-01-01-22. 
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APPENDIX A 


SYMBOL LIST 


A, a > 


V 

imaginary part of a complex 

B, b 

constants for complex 


quantity 

C, c 

► rational form and for 

WN 

a system constant that can be 

program functional form 


used in REDUCE 

D,d 

D 

denominator of G 

X 

signal or variable in a block 
diagram 

G 

transfer function for a block in 

X 

real part of a complex quantity 

i 

a block diagram 

YT 

y 

imaginary part of a complex 
quantity 

K 

a system constant that can be 

z 

a system constant that can be 

N 

used in REDUCE 
numerator of G 

z 

used in REDUCE 
a complex quantity 

S 

Laplace variable 

0) 

circular frequency, rad/sec 

T 

a system constant that can be 

Subscripts : 


used in REDUCE 

m 

index number 1,2,... 

u 

real part of a complex quantity 

n 

index number 1, 2, . . . 
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APPENDIX B 


USER'S MANUAL 

This appendix describes the use of both the FORMAC and FORTRAN programs pre- 
sented in this report. The description of the FORMAC program use is given first and 
the FORTRAN program second. The example used in this appendix is taken from the 
block diagram and system information given in figure 8 . 



(a) Block diagram. 


X 3 - GjX 2 * 0 

X 

no 

I 

X 

1— • 

X 

ON 

II 

0 

x 4 ' G 2 X 6 “ 0 

X 3 + X4 - X5 ■ 0 

Xi ■ - 0 


sinh(l) - sinh (d Vas 7 + bs ) 


cosh (U- cosh (dVas 2 + bs) 


VuT - Vas 2 + bs 


Solution required: x x /x 2 


(b) Algebraic information required. 

Kj * 763 359.38 

z 5 = 10.0 

Z x = 0. 75 

Z 6 =l-° 

Z 2 - 62. 4 

Tj = 0.010416 

Z 3 = 6.480x10 7 

a - 0.4100625x10 

Z4-32.2 

b = 0.0001 


d = 3. 33333 


(c) Numerical data required. 

Figure 8. - User's manual sample problem. 
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INPUT TO REDUCE 


Several points must be considered while preparing input to the FORMAC program. 

(1) All algebriac expressions must end with a $. 

(2) Blanks are ignored in the expressions and therefore may be used freely. 

(3) All 80 card columns of a card may be used for an expression. 

(4) An expression may take up to and including five cards. 

(5) The expressions of the block diagram may only be of the form 

±X(a)±X(b)±X(c) ... $ 


or 


±X(a)±X(b)* G(c)$ 


or linear combinations, as 

±X(a)±X(b)± . . . ±X(m)*G(j)±X(n)*G(k) ± . . . $ 

where a, b, c, etc. , are integer numbers. The number of terms in the first form must 
be greater than one. The =0 that makes the preceding expressions equations is implied 
and is not to be given. The total number of terms in all the block-diagram expressions 
may not be greater than 100. 

(6) The integer subscripts used in the X terms must start at one and increase se- 
quentially and must not be greater than 63. 

(7) The integer subscripts used in the G terms must start at one and increase se- 
quentially. Also, the largest subscript must not be greater than 20. 

(8) The numerator and denominator of the system information that is to replace the 
N's and D’s are input as separate expressions. These expressions are read sequen- 
tially; that is, the numerator for G(l), N(l), is read first, the denominator, D(l), is 
read second, the numerator for G(2) is read third, and so on. 

(9) The expressions of the system information may contain any of the following items: 

Item Description 

I V-L, i 

W omega, a; 

S iw 
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T(a),Z(a), 1 
K(a),WN(a) / 

EXPF(b) 

SINF(b) 

COSF(b) 

SINHF(b) 

COSHF(b) 

SQRTF(b) 

F(b) 

numbers 


variables 

complex exponential 
complex sine 
complex cosine 
complex hyperbolic sine 
complex hyperbolic cosine 
complex square root 
arbitrary function 

any number in FORTRAN notation, such as 1.0, 6.23E+2 


The subscripts a and b stand for integer numbers. These subscripts must begin at 
one and increase sequentially for each function and variable used. The maximum sub- 
script allowed is a = 10 (for variables) and b = 5 (for functions). 

Keeping the above points in mind, the input must be in the following order : 

(1) The first input card contains three integer numbers, (a) the number of variables 
(number of x’s), (b) the number of equations, and (c) the number of block-diagram in- 
puts. These numbers are punched on the card such that the last digit of the first number 
appears in card column 5, the second in card column 10, and the last in card column 15. 

(2) The second card contains X(a)$,X(b)$,X(c)$, . . ., where a,b, c, . . . are in- 
tegers that refer to the block-diagram inputs. 

(3) The third card contains X(a)$ and X(b)$, where a is the integer referring to 
the output variable of the transfer function desired and the integer b refers to the input 
variable of the transfer function. The ratio X(a)/X(b) is formed. 

(4) The block diagram expressions are punched on the next cards. One card may not 
contain more than one expression. 

(5) The word TRUE or FALSE is punched on the next card beginning in card column 
two. The word TRUE indicates that system information expressions are to follow. The 
word FALSE indicates that no system information follows, and items (6) and (7) are to be 
omitted from the input. 

(6) The system information expressions are punched on the next cards. One card 
may contain only one expression. 

(7) The word TRUE or FALSE is punched on the next card beginning in card 
column two. The word TRUE indicates that the program is to expand the transfer func- 
tion into real and imaginary parts. The word FALSE indicates that expansion is not 
desired. 
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Consecutive problems may be run by simply repeating the above steps. 

The following example of input to the FORMAC program is taken from figure 
8(a) and (b) : 


SOATA 

6 5 1 

X( 2 ) $ 

X ( 1 ) $ X ( 2 ) $ 

X( 2 )-X< 1 )-X (6)$ 

X(2)*G(1)-X(3)$ 

X(6)*G(2)-X(4)$ 

X ( 3 ) + X ( 4 ) —X ( 5 ) $ 

X( 5>*G( 3)-X( 1 ) $ 

TRUE 

— K ( 1 ) *Z ( 1 )**2*Z < 2 )*S**2*SINHF ( 1 )$ 
T ( 1 >*Z < 3)*Z <4)*SQRTF{ 1 )*COSHF ( 1 ) $ 
K ( 1 )+Z ( 5 ) *S$ 

1 $ 

1 $ 

Z ( 6 ) *S**2$ 

TRUE 


REDUCE OUTPUT 


Twelve errors are detected by the FORMAC program. Any of these errors will 
cause an error message to be printed and the termination of execution. Two warning 
messages may also be printed. The warning does not terminate execution. A list of the 
messages and their causes follow: 


Messages (errors) 

(1) MORE THAN 63 VARIABLES 

(2) MORE THAN 30 EQUATIONS 

(3) MORE THAN 5 BLOCK DIAGRAM 
INPUTS 

(4) MORE THAN 5 CARDS FOR A 
BLOCK DIAGRAM EXPRESSION 

(5) MORE THAN 5 CARDS FOR A 
SYSTEM INFORMATION EX- 
PRESSION 


Cause 

The number of variables exceeded the 
maximum allowed. 

The number of equations exceeded the 
maximum allowed. 

The number of block diagram inputs 
exceeded the maximum allowed. 

Block diagram expressions may not 
use more than 5 cards. 

System information expressions may 
not use more than 5 cards. 
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Messages (errors) 

(6) THE TOTAL NUMBER OF TERMS 
EXCEEDS 100 

(7) THE OUTPUT VARIABLE COULD 
NOT BE FOUND 


(8) THE INPUT VARIABLE COULD 
NOT BE FOUND 


(9) SYSTEM TOO LARGE 


(10) THE NUMBER OF TERMS IN- 
VOLVING X( ) IN EQUATION _ 
EXCEEDS 63 

(11) THE NUMBER OF TERMS IN- 
VOLVING X(_) IN ALL EQUA- 
TIONS EXCEEDS 63 

(12) SYSTEM CANNOT BE SOLVED 


Messages (warnings) 

(1) ONLY ONE TERM IN EQUA- 
TION 


(2) SOLUTION BELOW MAY NOT BE 
UNIQUE 


Cause 

The number of terms in all block dia- 
gram expressions exceeded the maxi- 
mum allowed. 

The output variable specified in the 
transfer function was not in the block 
diagram expressions. Check input 
equations . 

The input variable specified in the 
transfer function was not in the block 
diagram expressions. Check input 
equations . 

The internal array that contains the 
order in which the reduction is to pro- 
ceed was filled. 

During ordering the referenced varia- 
ble appeared in more than 63 reduc- 
tions in the equation referenced. 

During ordering the referenced varia- 
ble appeared in more than 63 reduc- 
tions . 

All variables were eliminated from the 
final equation. Check input equations. 

Cause 

The equation referenced contained only 
one variable during reduction. This 
may be a result of an error in the in- 
put equations. 

All block diagram equations were not 
used for the reduction. 


The printed output of the FORMAC program lists the following items: 

(1) The number of variables 

(2) The number of block-diagram equations 
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(3) The transfer function requested 

(4) The block-diagram equations after the variables that were specified as block- 

diagram inputs have been made zero 

(5) The variables that were specified as block-diagram inputs 

(6) The solution in the G-form 

(7) The solution in the N/D-form 

(8) The system information, if any 

(9) The solution in the S-form, if requested 

(10) The solution in the complex rational form, if requested 

Functional forms preceded by R or I indicate the real or imaginary parts asso- 
ciated with that function. Refer to equations (20) to (28). The punched output contains the 
solution in all the different forms. This output is punched in a form that is acceptable to 
the FORTRAN IV compiler. The S and complex rational forms are compatible with the 
FORTRAN programs given in this report. 

The examples of printed and punched output which follow were generated by REDUCE 
using the input listed in the Input to REDUCE section. 
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Printed Output from FORMAC Program 


THE MHBER OF VARIABLES IS 6 
THE NLMB ER OF EXPRESSIONS IS 5 


THE SULUT IUN IS FUR Milt 10 X(2)S 
THE ELQCK DIAGRAM EXPRESSIONS 

1 -XlllOIZl-XIflt 

2 G( 1 )* X I 2 )-X( 2 )$ 

3 G( 2 )*X<£ }-X( A )$ 

4 XI 3 J+XI4 )-X( f )$ 

5 G( 3)*X( 5 J~>( ] )S 


THE ELOCK DIAGRAM INPLIS 
X( 2 )$ 


N ( 

1 ) = 

-K( 1)*S**2.C*SINHF(I)*Z(1)**2.C*ZI2) i 

C{ 

1 )« 

CUSH F ( 1 I* SQR TF (1)*T(1)*Z!3)*Z(4)$ 

N ( 

2 ) = 

M 1)4<*Z( 5 )i 

01 

2 ) = 

1 .OS 

N ( 

3 ) = 

1 .cs 

C( 

3 ) = 

S**2.C*Z( £ )S 


THE DENCH IN AT OR 

G ( 2 ) *G( 3 ) + 1 .0 


THE NUMERATOR 

0tll*M2J*M3MCI2l*M 1I*N(3I 


THE 06 NCM IN AT OR 

D{1)*CI2)*DI 31*01 1I*M 2)*N(2) 


CCSHFI ll*CKC 1J4S*Z( 5 I )* SQR TF (1)* 1 1 l ) * Z I 3 )* 2 1 4 )- K (1 ) *S**2 ♦SINHFU 
I ♦Z C II **2 *Z I 2 ) 


Th£ “ISHMISFk. 1)*S*Z( S))*SQRTF( ll*mi*Z(3|*ZC4»4C0SHF (1)*S**2 * 

S CRT FI 1 > *T ( 1 »*Z( 3 )*Z( A )*ZI £ I 


THE NUMERATOR ( REAL I 


- I CCS H 1 ) * IS QRTC 1 J*Kl 1 )*T( 1 )*Z( 3 1* Z I 4 ) + K ( 1 J * RCOSH (1) ♦RSORT III *T (II 
*Z { 3 ) * 1 1 4 » 


KUJ*PSINHt 1 )*W*«2 *ZI1)**2 *Z I 2 ) 

(-ICCSHI 1 J*RSQRT I 1 1*1 ( 1 )*Z I 3 )*Z< 4 I* Z C 5)-I SQR Till *RC C SH 1 1 |»TI1J«2(3 
) *Z 1 4 ) *Z I 5 ) I *W 


THE NUMERATOR ( IM AG IN ARY I 


I CCS HI 1)*KI 1 I*RSQRTI 1 )*T( 1 )*Z13)*Z14) + I SQR T 1 1 > *K 1 l ) *RC C SH 1 1 ) *T 1 1 > * 
Z 1 3 ) *Z ( 4 ) 


ISINHI 1)*KI l)*W**2 *2111**2 * Z I 2 ) 


(-ICOSHI 1 )*ISQRT I 1 1*1 ( 1 l*ZI 3 l*ZI 4 l*ZI 5H-RC0SH (1) *RSQRT II ) *T 11 ) *Z (3 
J *Z I 4 I *Z I 5 ) ) *W 


THE DENOMINATOR (REAL) 


-IOCS HI 1 )*ISQRT( 1 )*KI 1 )*T( 1 )*Z I 3 I* Z I 4 ) ♦ K ( 1 ) *RCOSH ( 1 1 *RSQRTin*T(ll 
*Z I 3 ) *Z I 4 ) 


I ICOSM U + ISQRTI 1 )*T( 1 )*Z( 3 )*ZI 4 )* Z ( 6I-RC0SH ( 1) *RSQRT(I) *T(1 ) *Z(3) 
*Z ( 4 ) *Z I 6 ) )*W**2 


l-ICCSHI 1 )*RSQRT I 1 l*T ( 1 )*Z I 3 )* Z I 4 )* Z ( 5 )-I SORT (1 ) *RC CSH (1 I *T C 1 ) *Z (3 
)*Z(4)*Z( 5))*H 


THE DENOMINATOR (IMAGINARY) 


ICCSHI 1 ) *K( 1 )*RSQRT( 1 1*T( 1 I* Z ( 2 )*Z ( 4 ) ♦ I SQR T 1 1 ) *K (1 1 *RCOSH (1 ) *T 1 1 ) ♦ 
Z (3)*Z!4J 


(-ICCSH( 1 ) *RSQRT ( 1 )*T( 1 1*2 ( 3 )*Z ( 4 )* Z ( 6 )- 1 SQRT(1)*RC0SH (1 ) *T(1 ) *Z(3 

)*Z(4)*Z(6))*W**2 


l-ICCSHI I)*I$QRT( 1 J*T( 1 )*ZI 3 )*2 ( 4 )* Z I 5 )*RC OSH 1 1 ) *RSQRT ( 1 ) *T ( 1 ) *Z 13 
)*Z(4)*ZI 5) )*W 
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Listing of FORMAC Punched Output 


The G-form: 


C THE NUMERATOR 

NO = 

*G(1)#G(3)+G(2) *G I 3 ) 
C THE DENOMINATOR 

DO = 

*G ( 2 ) *G (31+1.0 


The N/D-form: 


C THE NUMERATOR 

N 1 = 

*D( 1 ) *N( 2 ) *N ( 3 ) +D( 2 ) *N ( 1 ) *N ( 3 ) 
c the denominator 

DI = 

*D( 1 )#D( 2 )*D(3 )+D( I )*N(2 )*N(3) 


The S-form: 


C THE NUMERATOR 

N2 = 

*COSHF ( 1 ) * ( K ( 1 ) +S#Z ( 5 ) ) *SQRTF ( I ) *T ( 1 ) *Z ( 3 ) *Z ( 4 ) -K ( 1 >*S**2 *SINHF( 1 
# ) *Z ( 1 ) **2 *Z<2) 

C THE DENOMINATOR 

D2 = 

*COSHF( 1 )*(K ( 1 )+S*Z ( 5 ) )*SQRTF < 1 ) *T < 1 ) *Z ( 3 ) *Z ( 4)+C0SHF ( 1 )*S**2 * 

#SQRTF ( 1 )*T( 1)*Z(3)*Z(4)*Z (6) 


The complex rational form: 


C THE NUMERATOR (REAL) 

A ( I ) = 

*- ICOSH ( 1 ) * I SORT < 1 )*K< 1 )*T(1)*Z (3) *Z ( 4)+K( 1 )*RCOSH( 1 )*RSQRT ( 1 )*T ( 1 ) 
#*Z ( 3 >*Z t A) 

A ( 2 ) = 

*K< 1 )*RSINH( 1 )*W**2 *Z ( 1 ) *#2 *Z(2) 

A ( 3 ) = 

* ( — I COSH ( 1 ) *RSQRT ( 1 ) #T ( 1 ) *Z ( 3 ) *Z ( A ) *Z ( 5 ) - 1 SOR T ( 1 )*RCOSH< 1 )*T( 1 )*Z <3 
*)*Z(4)*Z (5) )*W 
NA= 3 
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C THE NUMERATOR (IMAGINARY) 

B ( 1 ) = 

* ICOSH ( 1 )*K( 1 )*RSORT ( 1 )*T( I >*Z < 3 ) * Z ( 4 ) + I SOR T ( 1 ) *K ( 1 ) *RCOSH ( 1 ) *T ( 1 > * 
*Z(3)*Z(4) 

B ( 2 ) = 

*1 S INH ( 1 ) *K ( 1 ) *W**2 * Z ( 1 ) * * 2 *Z(2) 

B ( 3 ) = 

* ( - 1 COSH ( 1 ) * I SORT ( 1 ) * T ( 1 ) *Z ( 3 ) * Z ( 4 ) *Z ( 5 J+RCOSH ( 1 ) *RSO«T ( 1 ) *T ( 1 ) *Z ( 3 
*>#Z(4)*Z(5> >*W 
NB= 3 

C THE DENOMINATOR (REAL) 

C ( 1 ) = 

*— ICOSH ( 1 )*ISORT ( 1 )*K ( 1 )*T ( I )*Z ( 3) *Z (4)+K ( 1 )*RCOSH( 1 )*RSORT ( 1 ) ^ T ( 1 ) 
**Z(3)*Z<4) 

C( 2 ) = 

* ( ICOSH( 1 )*ISQRT ( 1 >*T< 1 >*Z (3 )*Z (4)*Z ( 6 )-RCOSH ( 1 ) *RSQRT ( 1 ) *T ( 1 ) *Z ( 3 ) 
**Z(4)*Z<6) )*W**2 
C( 3 ) = 

* ( - I C 0 S H ( 1 ) * R S 0 R T ( 1 ) * T ( 1 )*Z(3)*Z(4) > i : Z (5)-ISORT ( 1 )*RCOSH( 1 )*T ( I )*Z( 3 
*)*Z (4)*Z( 5) )*W 

MC= 3 

C THE DENOMINATOR (IMAGINARY) 

D ( 1 ) = 

*ICOSH( 1 ) *K ( 1 )*RSQRT ( 1 )*T ( 1 )*Z ( 3)*Z ( 4) + ISQRT( 1 )*K ( 1 )*RCOSH ( 1 )*T < 1 >* 
*Z ( 3)*Z( 4) 

D ( 2 ) = 

* (-ICOSH ( 1 )*RSQRT ( 1 )*T ( 1 >.*Z ( 3)*Z (4)*Z (6)-IS0RT( 1 )*RCOSH( 1 )*T ( 1 )*Z< 3 
*)*Z(4)*Z(6> )*W**2 

D ( 3 ) = 

*(-ICOSH( 1 )*ISQRT { 1 )*T< 1 )*Z ( 3)*Z (4)*Z ( 5)+RCOSH( 1 )*RSQRT( 1 )*T( 1 )*Z ( 3 
#)*Z(4)#Z(5) )*W 
ND= 3 


INPUT TO EVAL 

The user may evaluate the transfer function by using either the S-form or the com- 
plex rational form. If the S-form is to be evaluated, the S-form on punched cards is 
inserted between the two comment cards that contain asterisks in program COEFFS (see 
appendix C). The complex rational form requires the insertion of the punched cards be- 
tween the two comment cards that contain asterisks in program COEFFR. The numeri- 
cal results from both forms are identical. When the expression punched by EVAL ex- 
ceeds 19 continuation cards, the user must manually divide the expression into two or 
more parts. Each of these parts must contain less than 19 continuation cards to Comply 
with the FORTRAN IV compiler. 

If the user has made use of the arbitrary function in REDUCE, he must supply a 
function subprogram for his particular function. The numerical input consists of speci- 
fying values for the explicit variables used in the system information and values for the 
variables implied in the arguments of the allowed functions. A number corresponding to 
the maximum subscript employed in each variable and function used must be supplied. 
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The frequency range and increment desired must also be supplied. 

The first card of the EVAL input must contain $INPUT beginning in card column 2, 
and the last card must contain a $ in card column 2. The following input forms may ap- 
pear in any order but all must begin in card column 2: 

(1) NEXP=number, EXPA=list, EXPB=list, EXPC=list, EXPD=list, 

(2) NSIN=number, SINA=list, SINB=list, SINC=list, SIND=list, 

(3) NCOS=number, COSA=list, COSB=list, COSC=list, COSD=list, 

(4) NSINH =number , SINHA=list, SINHB=list, SINHC=list, SINHD=list, 

(5) NCOSH=number, COSHA=list, COSHB=list, COSHC=list, COSHD=list, 

(6) NSQRT=number, SQRTA=list, SQRTB=list, SQRTC=list, 

(7) NK=number, K=list, 

(8) NT=number, T=list, 

(9) NZ=number, Z=list, 

(10) NWN=number, WN=list, 

(11) FSTART=initial frequency, FEND=final frequency, DELTA F=increment, 

Any of these forms may use more than one card as long as each card ends with a 
comma. The suffixes A, B, C, and D in the function names represent the variables 
implied in the function argument. (Note: no D inSQRT.) For example, EXPA, EXPB, 
EXPC, and EXPD represent the variables a, b, c, and d in the function 


exp \d Vas 2 + bs + c 

(See eqs. (15) to (20) for additional forms.) The term "number, ” as used in the pre- 
ceding forms, means an integer number which has the value of the maximum subscript 
used for that function or variable. For the functions, this number must not be greater 
than 5; for variables, not greater than 10. The term "list" means a list of numbers that 
the variables are to take on. The elements of this list must be separated by commas and 
must be in the proper order. For example, NK=2, K=l, 6. 25E-2 would set K(l)=l. 0 and 
K(2)=0.0625. For consecutive runs (runs where the transfer function does not change) 
only the changed numbers need to be given on the input. For example, 

SINPUT 

K ( 2 ) =1 .25E-2 » 

$ 
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In the following example of input to EVAL, it is understood that either the S-form of 
the complex rational form of the transfer function (given in the FORMA C Sample Output) 
from REDUCE is to be inserted in the proper FORTRAN subroutine. These numerical 
data are taken from figure 8(c). 


$DAT A 
SINPUT 

NSINH=1, SINHA=4.100625E— 6,SINHB=1.0E-4,SINHC=0,SINHD=3.33333, 
NC0SH=1 ,C0SHA=4. 10062 5E-6,COSHB=l .0E-4,C0SHC=0 t C0SHD=3. 33333, 
NS0RT=1 ,SQRTA=4. 10062 5 E-6, SQRTB=1 .OE-4, SQRTC=0, 

NK=1 ,K=.76335938E6, 

NT = 1 y T=1.0416E-2, 

NZ =6,Z =0.750, 62. 4,6. 48 E7, 32. 2,10,1, 

FSTART=1 , FEND=50,OELTAF=1 t 
$ 


EVAL OUTPUT 

The output printed by the FORTRAN program lists the following items: 

(1) All numerical input 

(2) For each frequency requested 

(a) The frequency 

(b) The real and imaginary parts of the numerator and denominator of the trans- 

fer function normalized so that either the real or imaginary part of the 
numerator is one 

(c) The real and imaginary parts of the transfer function 

(d) The complex absolute value of the transfer function 

(e) The complex amplitude (phase angle) of the transfer function 

An example of the EVAL output follows. This output was obtained by using the input 
given in the previous section. 
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FORTRAN EVALUATION OF FORMAC OUTPUT 


I 

SIN LA 

SINHB 

SI NHC 

SINHD 

1 

0.41CC6E-05 

C. 100C0E-03 

0 

3.33333 

I 

COSLA 

CQSHB 

CCSHC 

COS HO 

1 

0.41CC6E-C5 

C.10000E-03 

0 

3.33333 

I 

SQR1A 

SQRTB 

SCRTC 


1 

0.4 1CC6E-05 

C.10000E-03 

0 


I 

K 




1 

0.76236E 06 




I 

T 




1 

0.10* 16E-0 1 




I 

l 




1 

C.75CCC 




2 

62.4CC0 




3 

0.648C0E C 8 




4 

22.2CCO 




5 

10.CCC0 




6 

1 .CCCOO 





FSTART 

= 1 .COCCO 

FEND= . 

5C.0C00 

DELT AF= 

1 .00000 





p 

A 

B 

C 

0 

REAL 

I MAG 

AB 

1 

l.OOOCC 

1.00000 

1.29988 

0.99974 

1.29954 

1.00026 

-0. 52105E-06 

1.00026 

2 

2 .COOCC 

l.CCOOO 

1. 66699 

0.99893 

1.66523 

1 ,0C 106 

- C • 41 5 7 IE- 05 

1. 00106 

3 

3.00CCC 

1.00000 

2. 09351 

0.99759 

2.08853 

1 .CC 23 S 

—0.141 19E- 04 

1. 0023 9 

4 

4. COOCC 

l.COOOO 

2. 57234 

0.99567 

2.56145 

1 .00 A 26 

- 0 . 33830E- 04 

1.00426 

5 

5.00CCC 

l.COOOO 

3.10066 

0.99314 

3.080 10 

1 .0C665 

- 0 .66 9 4 IE- 04 

L. 00669 

6 

6.000CC 

l.COOOC 

3. 68043 

0.98997 

3.64520 

1.CC970 

-C.il 756E-03 

1. 00970 

7 

7 .000 CO 

l.CCOOO 

4.31822 

0.98608 

4.26174 

1.01229 

-C. 19034E-03 

1.01329 

a 

8.00CCC 

l.OOCOO 

5. C2519 

0.98138 

4.93900 

1.01751 

- C • 29062E- 03 

1. 01751 

9 

9. 00 C CO 

l.COOOO 

5. 81798 

0.97575 

5.69107 

1.02237 

-0.42466E-03 

1. 02237 

10 

1 C .occc 

1. COOOO 

6. 72 03 9 

0.96902 

6.53839 

1.C2792 

-0 . 59990E- 03 

1. 02792 

11 

il.OCCC 

1.00000 

7. 766 57 

0.96093 

7.51047 

1.02420 

-C.82526E-03 

1. 0342 0 

12 

12.00CC 

l.COOOO 

9. 00666 

0.95114 

8.65074 

1.04127 

-0.1111 6E -02 

1.04127 

13 

13. OCCC 

1. CCCOO 

10.5169 

0.93907 

10 .0253 

1.C4917 

-0.1472 IE- 02 

1.04917 

14 

14.00CC 

1.00000 

12.4180 

0.92384 

11.7398 

1.05759 

- G . 19232E- 02 

1. 05800 

15 

15. OCCC 

l.COOOO 

14.9189 

0.90396 

13.9734 

1.C6782 

-C.24854E-02 

1. 06782 

16 

16 .00 CC 

l.COOOO 

1 8.3976 

0.87664 

17.0570 

1 .C 7 876 

-0.31 847E- 02 

1. 07876 

17 

17 .00 CC 

l.COOOO 

23.6435 

0.83610 

21.6759 

1 .05053 

-C.40535E-02 

1.09094 

18 

18.00CC 

l.COOOO 

32.5991 

0.76814 

29.5185 

1 • 1C 445 

-C.51356E-C2 

1.10450 

19 

19. OCCC 

l.COOOO 

51.6918 

0.62572 

46.1724 

1.11562 

-C.64851E-02 

1. 11964 

20 

2C.00CC 

l.COOOO 

122.467 

0.10475 

107.753 

1.12656 

-C.01755E-O2 

1.13659 

21 

21.C0CC 

l.COOOO 

-333. 826 

3.44122 

-288.852 

1.15558 

-0.103C5E-01 

1.15562 

22 

22. OCCC 

l.COOOO 

-7 C. 0887 

1.50742 

-59.5307 

1.177C2 

- C.130C6E-01 

1.17710 

23 

23. OCCC 

l.COOOO 

-38.7821 

1.27455 

-32.2645 

1.2C136 

- 0 . 16464E- 01 

1. 20147 

24 

24.C0CC 

l.COOOO 

-2 6.53 0 4 

1 .18090 

-21.5644 

1.22514 

-0.20937E-01 

1.22932 

25 

25.00CC 

1. COOOO 

-1 9.9393 

1.12841 

-15.7068 

1.26112 

- 0 . 26 7 99E- 01 

1. 26140 

26 

26.00CO 

l.COOOO 

-15. 7847 

1.09352 

-12.1290 

1 .25826 

— C . 34603E- 01 

1.29875 

27 

27 .CCCC 

l.OOCOO 

-12.9000 

1.06766 

-9 .57693 

1.24195 

-C.45187E-01 

1.34271 

28 

28.C0CC 

l.COOOO 

-1 C. 7606 

1.04701 

-7.67463 

1.35254 

— 0 , 59868E- 01 

1.39522 

29 

29.O0C0 

l.COOOO 

-9. 09548 

1.02959 

-6.18662 

1.45674 

-0.80753E-01 

1. 45898 

30 

30.00CC 

l.CCOOO 

-7. 75048 

1.01426 

-4.97902 

1.52285 

-0.11162 

1.53794 

31 

31. OCCC 

l.COOOO 

-6. 63146 

1.00034 

-3.97007 

1 .6 3C 3 2 

-C. 15 891 

1 . 63805 

32 

32.00CC 

1.00000 

-5.67750 

0.98738 

-3.10700 

1.75264 

-0.23512 

1. 76834 

33 

33. OCCC 

l.COOOO 

-4. 84777 

0.97509 

-2.35419 

1 .SC 78 5 

-0.36544 

1. 94253 

34 

34.00CO 

1.00000 

-4.11324 

0.96325 

-1.68675 

2.C5418 

-0.60306 

2.17928 

35 

35. OCCC 

1.00000 

-3.45320 

0.95171 

-1.08672 

2.25445 

- 1.05416 

2.48073 

36 

36.00CC 

l.COOOO 

-2. 85222 

0.94036 

-0 . 540 8 1 

2. 1C554 

- 1.81967 

2.78623 

37 

37.C0CC 

1.00000 

-2.29859 

0.92911 

—0 .389756-01 

1.178C1 

-2.42455 

2. 69558 

38 

38.C0CC 

1.00000 

-1.78321 

0.91790 

0 .42655 

C. 15352 

-2.01405 

2. 01989 

39 

39. OCCC 

l.COOOO 

-1.29891 

0.90667 

0 .06185 

-0 .12555 

- 1.30335 

1.31043 

40 

40.00CC 

1. 00000 

-0. 83991 

0.89538 

1.27101 

-0 .71441 E-01 

-0.03658 

0.83962 

41 

41 .OCCC 

l.COOOC 

-0.40149 

0.88399 

1.660 36 

G. 61437 E- Cl 

-C .56957 

0. 57287 

42 

42.00CC 

l.CCOOO 

0. 203 03 E-01 

0.87248 

2.03075 

0.18704 

-C. 41207 

0.45253 

43 

43.00CC 

l.COOOO 

0. 42878 

0.86081 

2.38568 

0.25285 

-C. 31350 

0.42900 

44 

44.00CO 

l.COOOO 

0. 82680 

0.84897 

2.72744 

C .38041 

-C. 24823 

0.45424 

45 

45 .00 CC 

l.COOOC 

1.21687 

0.83692 

3.05799 

C. 452*7 

-C. 20291 

0.49679 

46 

46.00CC 

1.00000 

1.60122 

0.82466 

3.37901 

0.5154C 

-0.17016 

0. 54276 

47 

47 .OOCC 

1. CCCOO 

1. 981 84 

0.81216 

3.6920 2 

C .56805 

-0 . 14572 

0. 5 872 L 

46 

48.00CC 

l.OOCOO 

2.36061 

0.79941 

3.99836 

C .61578 

-C. 12695 

0.62874 

49 

49 .OOCC 

1. COCOO 

2. 73926 

0. 78639 

4.29923 

C. 65769 

-C. 11 230 

0.66721 

50 

50. OCCC 

1. CCCOO 

3.11947 

0.77307 

4.59578 

C. 65565 

-C. 10057 0. 702 92 

REC = 00000 F I L = 


ANGLE 

-0.29846E-04 

-0.23793E-03 

-0.80704 E-03 

-0.19301E-02 

-0.38099 E-02 

-0,667 12 E-02 

-0.10763 E-0 1 

-0.163656-01 

-0.23799 E-01 

-0.33438E-01 

-0.457206-01 

-0.61 164E— 01 

-0.80391E-01 

-0.10415 

-0.13336 

-0.16915 

-0.21291 

-0.26641 

-0.33186 

-0.41214 

-0.51092 

-0.63310 

-0.78514 

-0.97585 

-1.21735 

-1.52672 

-1.92856 

-2.45929 

-3.17447 

-4.16203 

-5.56699 

-7.64072 

-10.8434 

-16.0647 

-25.0604 

-40.7754 

-64.0865 

-85.6412 

-95.9567 

-94.8811 

-83.8435 

-65.5869 

-46.9507 

-33.1263 

-24.1065 

-18.2707 

-14.3684 

-11.6521 

-9.68961 

-8.22560 
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APPENDIX C 


DESCRIPTIONS, LISTINGS, AND FLOW CHARTS FOR COMPUTER PROGRAMS 

FORMAC Program 

REDUCE (main program) . - This FORMAC program is the body of the FORMAC 
programs. It reads the input, reduces the equations, performs the substitutions, ex- 
pands and separates the real and imaginary parts, and calls other subroutines for ini- 
tialization and output. 

Subrout ine SOUT . - This FORMAC subroutine writes and punches the symbolic 
transfer function. To do this, SOUT calls the subroutine that removes the $ and floating 
point powers from the symbolic expressions. 

Subroutine INIT, subrout ine SORT . - The purpose of subroutine INIT is to determine 
the order in which the block diagram equations are to be reduced. For a detailed de- 
scription of the method used, see appendix D. 

Subroutine SORT sorts an array of integer numbers into an array of numbers that in- 
crease in magnitude . Subroutine SORT is used by subroutine INIT . 

Subroutine COS, SIN, THNH, ATAN, ALOG, DUMP, PDUMP, FMCDIF, EX DDMP , 
FMCDMP . - This is a dummy subroutine that is not executed. Its purpose is to provide 
entry point names that are the same as the subprogram's that are referenced by FORMAC, 
but not executed by the programs in this report. This causes the programs that are not 
needed not to be loaded. 

Subroutine REMOVE, subroutine DELETE, function NBR, function FC. - This set of 
subprograms removes the $ and floating point powers from the BCD representation of 
FORMAC expressions. They are called just prior to the punching of any expression. 


FORTRAN Program 

EVAL (m ain program) . - Program EVAL reads and writes all the numerical input 
and output. It also calls the subprogram that evaluates the symbolic transfer function. 

Function EXPFX, SINFX, COSFX, SINHFX, COSHFX, SQRTFX. - This function sub- 
program evaluates the allowed system functions. Different entry point names are used 
for each function. 
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Subroutine COEFF. - Two subroutines with the name COEFF are supplied. The deck 
names are used to differentiate between them. The subroutine with the deck name 
COEFFS is used for transfer functions in the S-form. Deck COEFFR is used for the 
complex rational form of the transfer function. In both cases, the symbolic transfer 
function is to be inserted between the two comment cards that contain asterisks . 

The rest of this appendix gives listings and flow charts (figs. 9 to 14) for the pro- 
grams. Flow charts are only given for the programs in which the flow warrants a chart. 
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Figure 9. - Flow Chart of REDUCE. 
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Computer Listings 
REDUC 


$IBFMC REDUCE 
SYMARG 

ATOMIC X(63>,G(20),N(20),D(20),I,T(10),Z(10),K(10),WN(10),W,S, 
*EXPF<5) ,REXP<5) , l£XP<5> ,SINF(5) , RSI N < 5 > , I S I N( 5 ) ,COSF < 5 ) ,RCOS<5) * 
*ICOS( 5) , SINHF( 5 ) ,RS I NH < 5 ) , ISINHI 5 ) ,COSHF ( 5) ,RCOSH( 5) , I COSH! 5) , 
*SQRTF(5 ) tRSORT ( 5 ) , I SORT ( 5 ) , F 1 5 ) , RF ( 5 ) , I F ( 5 ) 

DIMENSION I N ( 60 ) ,OUT < 1 7 ) , NEQ ( 1 50 > , EQ ( 30 ) , INPUT) 5 > , SN ( 20 ) , SD ( 20 ) 
EQUIVALENCE t I N ( 1 ) ,OUT< 1 ) ,NEQ( 1 ) ) 

LOGICAL LOCATE, TEST, REPLCE,EXPND 

REAL I ,K,N, IS, IT, I EXP, I SIN, I COS, I S I NH , I COSH , NUM, I NPUT , I SORT , I F 
CMPX PARAM ( I**L»I >, ( I**(L+1) ,-l) ,< I**(L+2),-I ) , < I**(L+3) ,1 ) 

FNCT PARAM ( EXP F < J ) , R E X P < J ) + 1 * I EXP ( J ) ) , 

* ( SINF) J ) ,RS IN< J )+I*ISIN( J ) ) , 

* (COSF(J) ,RCOS( J)+I*ICOS(J ) ) , 

* < SINHF( J ) ,RSINH( J) + I*ISINH< J ) ) , 

* ( COSHF ( J ) , RCOSH t J ) + I * ICOSH ( J ) ) » 

* ( SQRTF ( J ) ,RSQRT< J ) + I*I SQRT( J ) ) , 

* <F< J) ,RF(J)+I*IF( J) ) 

DATA DR /6H $00000/ 

C READ ALGEBRAIC EQUATIONS 

£*****$£**#* #**#*#*>!<***********$*******$**£***’('****>!'********=!'*>!<:!<******** 

98 READ (5,101) M,M1,IK 
IFIM.GT.63) GO TO 401 
IFIM1.GT.30) GO TO 402 
I F ( IK.GT.5) GO TO 403 
JJ = 0 

READ (5,102) ( IN( J ) , J = 1 , 12 ) 

DO 96 J=1,IK 

LET INPUT(J) = ALGCON IN(1),JJ 

96 CONTINUE 
JJ = 0 

READ (5,102 ) ( I N ( J ) , J = l, 12 ) 

LET NUM = ALGCON IN(1),JJ 

LET DEN = ALGCON IN(1),JJ 

DO 100 L = 1 , M 1 

JJ = 0 

DO 99 J=1 ,5 

K 1 = 12 * ( J— 1 ) + 1 

K 2 = Kl + 12-1 

READ (5,102) ( IN(KK > ,KK=K1,K2) 

DO 99 KK = 1 ,72 

99 I F ( FC ( IN(K1 ) ,KK ) . EQ.DR ) GO TO 97 
WRITE (6,230) 

STOP 

97 LET EO ( L ) = ALGCON IN(1),JJ 
LET EQ ( L ) = EXPAND EQ ( L ) 

100 CONTINUE 

INITIALIZE AND SET UP EQUATIONS FOR SOLUTION 

*#*#**£ St*^^^*^*** #*#**«#*## *c J*#*# 

FIND THE NUMBER OF FUNCTIONS USED 

DO 29 L= 1 , 20 
DO 27 J = 1 , M 1 

LET LOCATE = FIND EQ ( J ) , APP , ALL , ( G ( L ) ) 

27 I F ( LOCATE ) GO TO 28 
GO TO 30 

28 MG = L 


45 


no non non non 


29 CONTINUE 
C 

C REPLACE THE INPUTS BY ZERO IF NOT DEN OR NUM 

C 

30 DO 22 L=1,M1 
DO 21 J=1 » I K 

LET TEST = MATCH I D, NUM , I NPUT < J > 

I F ( TEST ) GO TO 21 

LET TEST = MATCH I D, DEN , INPUT < J ) 

I F ( TEST ) GO TO 21 

LET EO ( L ) = SUBST EO < L ) , ( I NPUT { J ) , 0 ) 

21 CONTINUE 

22 CONTINUE 

£444444 44444444444444444444 44 4*#4444*#44#444444*444444444444444444444444 

C WRITE ALGEBRAIC EQUATIONS 

£44444444444444444444444444444444444**4*44444444444444444 444444444444444 

WRITE (6,201) M , Ml 

Q = 0.0 

LET Q = BCDCON NUM, OUT, 17 
WRITE (6,202) OUT ( 2 ) 

Q = 0.0 

LET 0 = BCDCON DEN, OUT, 17 
WRITE (6,203) 0UT(2) 

WRITE (6,204) 

DO 301 L = 1 , M 1 

Q = 0.0 

51 LET Q = BCDCON EQ(L),0UT,10 

WRITE (6,208) L, (OUT(J) ,J=2,10) 

IF(Q.NE.O.O) GO TO 51 
301 CONTINUE 

24 WRITE (6,213) 

DO 25 L=1,IK 

0 = 0.0 

LET Q = BCDCON I NPUT ( L ) ,OUT , 17 

25 WRITE (6,206) 0UT(2) 

£44444444444444444444444444444444444444444444444444444444444444444444*44 

C REDUCE THE EQUATIONS 

26 CALL INIT (NEQ,X»M,EQ,M1,NUM,0EN) 

GET A VARIABLE TO ELIMINATE 

11=0 
611= II+l 
J = NEQ ( I I ) 

IF(J.EQ.O) GO TO 12 

GET AN EQUATION CONTAINING THE VARIABLE 

6 II = II+l 
L = NEQ ( I I ) 

SOLVE THE EQUATION 

8 LET RR = EQ ( L ) 

ERASE EO ( L ) 

LET IT = COEFF RR,X(J) 

LET RR = EXPAND RR-IT*X(J) 

ELIMINATE THE VARIABLE IN THE REMAINING EQUATIONS 


non 


c 

9 II = II+l 
L = NEQ< I I ) 

IF(L.EQ.O) GO TO 10 
LET IS = COEFF E0(L),X(J) 

LET £Q( L ) = EXPAND EO ( L ) -I S*X ( J ) 

LET EO ( L ) = EXPAND I S*RR- I T*EQ ( L ) 

GO TO 9 

10 ERASE IS, IT 
GO TO 4 

REMOVE FUNCTIONS RAISED TO A POWER IF POSSIBLE 

12 II = II+l 
L = NEQ ( I I ) 

LET R = EO ( L ) 

ERASE EO ( L ) 

DO 13 L=1 , MG 
J = 0 

LET IS = COEFF R ,G ( L ) ** J , V 1 , V2 
LET TEST = MATCH ID, IS, 0.0 
ERASE IS 

I F ( .NOT. TEST) GO TO 13 
LET R = EXPAND R*G ( L ) ** ( -V 1 ) 

13 CONTINUE 

C * * ****** * ******** Jjc# ## *«#*#** ajc Jit*#*## ***#:{(*:*#* if $ * * $***#*# jft * 

C SOLVE FOR NUMERATOR AND DENOMINATOR 

LET RD = COEFF R , NUM 
LET RN = COEFF R , DEN 
LET RD = EXPAND -1*RD 

C $*#***>)c**#$:j( **#**))■#* ** #****#*# A##***####*^#*#:*;*********#*****###* 

C WRITE NUMERATOR AND DENOMINATOR 

WRITE (6,205) 

PUNCH 221 

CALL SOUT (RN,W,2,2HN0) 

WRITE (6,207) 

PUNCH 222 

CALL SOUT (RD,W,2,2HD0) 

ERASE RN , RD 

C SUBSTITUTE • N ( L ) / D ( L ) • FOR 'G(L)' AND REMOVE D(L)**(-1) 

DO 14 L = 1 , MG 

LET LOCATE = FIND R , APP , AL L , ( G ( L ) ) 

IF( .NOT. LOCATE) GO TO 14 
J = 0 

LET IS = COEFF R ,G ( L ) **J , V 1 , V2 
ERASE IS 

LET R = SUBST R,(G(L),N(L)/0(L)) 

LET R = EXPAND R*D(L)**V2 

14 CONTINUE 

C SOLVE FOR NUMERATOR -AND DENOMINATOR 

C * * * # * # * * * $ * # * # * 4 * ## # £ * £ * # $ $ 4 * if if if if if if if >!= 4 if if if if if if if >1= $ if if if if * if if * * * if if jjt $ * if if # ** ♦ # * # 
LET RD = COEFF R,NUM 

LET RN = COEFF R , DEN 

LET RD = EXPAND -1*RD 

ERASE R 
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C*********************************************************************** 

C WRITE NUMERATOR AND DENOMINATOR 

C**#**********:********* **************** ********************************* 

WRITE (6,205) 

PUNCH 221 

CALL SOUT (RN,W,2,2HN1 ) 

WRITE (6,207) 

PUNCH 222 

CALL SOUT (RD,W,2,2HD1) 

£*£********************************************************************* 
C READ AND WRITE THE SYSTEM INFORMATION 

£#!!<*#** ************************** *************************************** 
READ (5,103) REPLCE 
I F ( .NOT.REPLCE) GO TO 98 
II = 0 
MG2 = 2*MG 
DO 61 L= 1 , MG2 
LL = MOD ( L ,2 ) + 1 
JJ = 0 
DO 58 J=1 ,5 
K 1 = 1 2* ( J-l ) +1 
K2 = Kl+12-1 

READ (5,102) ( IN(KK) ,KK=K1,K2) 

DO 58 KK=1 ,72 

58 IF ( FC( IN(K1 ) ,KK ) .EO.DR ) GO TO (60,59),LL 
WRITE (6,231) 

STOP 

59 II = 11+1 

LET SN(I1) = ALGCON IN(1),JJ 
GO TO 61 

60 LET SO ( II) = ALGCON IN(1),JJ 

61 CONTINUE 

0 = 0.0 

DO 34 L= 1 , MG 
WRITE (6,214) L 

32 LET 0 = BCDCON SN(L),0UT,9 

WRITE (6,215) ( OUT ( J ) , J=2 , 9 ) 

IF(Q.NE.O.O) GO TO 32 
WRITE (6,216) L 

33 LET 0 = BCDCON SD(L),0UT,9 

WRITE (6,215) ( OUT ( J ) » J = 2 , 9 > 

IF(O.NE.O.O) GO TO 33 

34 CONTINUE 

£**#*#>!:)!!!!!:(£:{=#******************* ********************** ************ ******* 
SUBSTITUTE THE SYSTEM INFORMATION FOR * N ( L ) • AND • D ( L ) • 
*********************************************************************** 
DO 35 L=1 , MG 

LET RN = SUBST RN , ( N ( L ) , SN ( L ) ) , ( D ( L ) , SD ( L ) ) 

LET RD = SUBST RD, ( N ( L ) , SN ( L ) ) , ( D ( L ) , SD ( L ) ) 

ERASE SN ( L ) , SD ( L ) 

35 CONTINUE 

C************************* ************************** ******************** 

C WRITE NUMERATOR AND DENOMINATOR 

C**#*****##* **************************************************** ******** 

WRITE (6,205) 

PUNCH 221 

CALL SOUT (RN»W,2,2HN2 ) 

WRITE (6,207) 

PUNCH 222 

CALL SOUT (RD,W,2,2HD2 ) 



C READ EXPND 

C *£## * £ #£«*# £>!>* & * >1= if if if if if if if if if if if if *#:(:# s* ## # « >je **£ ajc ajt $ 

READ (5,103) EXPND 

IF( .MOT. EXPND) GO TO 98 

C * # * if * * if * # # * * if if if # if if # * *e if * * =!< # if if * * * if * ^ $ s|t if if sjc £ if * * * * if * * ❖ * * # * * * if * # * * « * # « « sit if if if if 

C SUBSTITUTE REAL FORMS FOR COMPLEX FORMS AND EXPAND 

C * * * * #*** * ** ** ** * if * * * * ****£# * * =!c ^ if if £ ** =!c sjt if if if if if ** if if * if £ =jc j)t # if *ififififif *$** **#**# * 

DO 38 J = 1 » 5 

LET RN = SUBST RN , FNCT 
LET RD = SUBST RD , FNCT 
38 CONTINUE 

LET RN = SUBST RN,(S,W*I) 

LET RD = SUBST RD,(S,W*I) 

LET RN = EXPAND RN 
J = 0 

LET IS = COEFF RN, I **J ,V1 , V2 
ERASE IS 
K 1 = V2 

IFIK1.LT. 2) GO TO AO 
DO 36 L=1,K1,4 
LET RN = SUBST RN,CMPX 
36 CONTINUE 

C SEPERATE INTO REAL AND IMAGINARY PARTS AND WRITE 

40 LET B 1 = COEFF RIM, I 

LET A 1 = EXPAND RN-I#B1 
ERASE RN 
WRITE (6,217) 

PUNCH 223 

CALL S0UT(A1,W,1,1HA) 

ERASE A 1 
WRITE (6,218) 

PUNCH 224 

CALL S0UT(B1,W,1,1HB) 

ERASE B 1 

LET RD = EXPAND RD 
J = 0 

LET IS = COEFF RD , I ** J , V 1 , V2 
ERASE IS 
K 1 = V2 

IFIK1.LT. 2) GO TO 42 

DO 41 L = 1 »K1 ,4 

LET RD = SUBST RD,CMPX 

41 CONTINUE 

42 LET D1 = COEFF RD, I 

LET Cl = EXPAND RD-I*D1 
ERASE RD 
WRITE (6,219) 

PUNCH 225 

CALL SOUT (C1,W,1,1HC) 

ERASE Cl 
WRITE (6,220) 

PUNCH 226 

CALL SOUT (D1,W,1,1HD) 

ERASE D1 
GO TO 98 

C ** if ^ if if if if if ^ if 4 if if if if if if if if if if if if if if if * $ * * if * if if * if * * if ## if ifififif $ # # * ***« * i if if Jjc if if if if if * if * * if * 

401 WRITE (6,227) 




GO TO 500 

402 WRITE (6,228) 

GO TO 500 

403 WRITE (6,229) 

500 STOP 

101 FORMAT (315) 

102 FORMAT (13A6,A2) 

103 FORMAT (L5) 

201 FORMAT (1H1,//,27H THE NUMBER OF VARIABLES IS,I5,/,29H THE NUMBER 
*0F EXPRESSIONS IS, 15) 

202 FORMAT ( 21HKTHE SOLUTION IS FOR ,A6,3HT0 ) 

203 FORMAT (1H+,29X,A6) 

204 FORMAT ( 30HKTHE BLOCK DIAGRAM EXPRESSIONS,//) 

205 FORMAT ( 1 HK , 55X , 1 3HTHE NUMERATOR) 

206 FORMAT (5X,9A6) 

207 FORMAT ( 1HK , 55X , 1 5HTHE DENOMINATOR) 

208 FORMAT ( 3X , 1 2 , 2 X , 9A6 ) 

213 FORMAT ( 25HKTHE BLOCK DIAGRAM INPUTS,//) 

214 FORMAT ( 3H N(,I2,2H)=) 

215 FORMAT ( 1H+ , 7X , 8A6 , / / ) 

216 FORMAT ( 3H D(,I2,2H)=) 

217 FORMAT ( 1HK ,55X ,20HTHE NUMERATOR (REAL)) 

218 FORMAT ( 1HK,55X,25HTHE NUMERATOR (IMAGINARY)) 

219 FORMAT ( 1HK ,55X ,22HTHE DENOMINATOR (REAL)) 

220 FORMAT ( 1 HK , 55X , 27HTHE DENOMINATOR (IMAGINARY)) 

221 FORMAT ( 1HC ,5X » 13HTHE NUMERATOR) 

222 FORMAT ( 1 HC , 5X , 1 5HTHE DENOMINATOR) 

223 FORMAT ( 1 HC , 5X , 20HTHE NUMERATOR (REAL)) 

224 FORMAT ( 1 HC , 5X , 2 5HTHE NUMERATOR (IMAGINARY)) 

225 FORMAT ( 1 HC , 5X , 22HTHE DENOMINATOR (READ) 

226 FORMAT ( 1 HC , 5X , 27HTHE DENOMINATOR (IMAGINARY)) 

227 FORMAT ( 33H TERROR** MORE THAN 63 VARIABLES) 

228 FORMAT ( 33H EEERRORE* MORE THAN 30 EQUATIONS) 

229 FORMAT (43H EEERROR** MORE THAN 5 BLOCK DIAGRAM INPUTS) 

230 FORMAT (59H **ERROR** MORE THAN 5 CARDS FOR A BLOCK DIAGRAM EXPRES 

*SION) 

231 FORMAT (64H E*ERROR** MORE THAN 5 CARDS FOR A SYSTEM INFORMATION E 
EXPRESSION) 

END 
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SOUT 


$ I 8FMC SPLIT 

SUBROUTINE SOUT < X , W, I , N AME ) 
DIMENSION OUT (12) 

EQUIVALENCE ( OUT ( 1 ) , MOUT ) 

LOGICAL TEST 

SYMARG X,W 

ATOMIC W 

GO TO (4,5 ) , I 

4 A = 0 
J = 0 

1 LET Y = COEFF X,W**A,V1,V2 
LET Y = Y*W**A 

A = V2 

LET TEST = MATCH ID, Y, 0.0 
IF(TEST.AND.A.EQ.O.O) GO TO 7 
I F ( TEST ) GO TO I 
WRITE (6,203) 

J = J + l 

PUNCH 204, NAME, J 

2 LET 0 = BCDCON Y, OUT, 12 
CALL R EMOVE ( OUT ) 

3 WRITE (6,201) ( OUT ( L ) , L = 2 , 12 ) 

PUNCH 202 , ( OUT (L)»L = 2,12) 
IF(O.NE.O.O) GO TO 2 

ERASE Y 

IF(A.NE.O.O) GO TO 1 
7 PUNCH 205, NAME, J 
RETURN 

5 PUNCH 206, NAME 

6 LET Q = BCDCON X, OUT, 12 
CALL REMOVE ( OUT ) 

WRITE (6,201) ( OUT ( L ) , L = 2 , 12 ) 

PUNCH 202, (OUT (L) ,L=2 ,12) 
IF(Q.NE.O.O) GO TO 6 
RETURN 


201 

FORMAT 

(61X,11A6) 

202 

FORMAT 

(5X,1H*,11A6) 

203 

FORMAT 

( 1HK ) 

204 

FORMAT 

(6X,A1,1H( , I 2 ,2H ) = 

205 

FORMAT 

( 6 X , 1 HN ,A1,1H=, 12) 

206 

FORMAT 

END 

<6X,A2,1H=) 
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INIT 


$ I BFMC INIT. 

SUBROUTINE IN I T < NEO, X , M , EQ , Ml ,NUM , DEN ) 

DIMENSION NEO (150) ,EQ<30) ,N0( 100) 

LOGICAL LOCATE, 0UT1,0UT2 
REAL NUM 

INTEGER DN,S,SS,SSS,S2,S3,S4,S5,S6,EX4,EX5,EX6,EX 
DATA S2, S3 ,S4 , S5 , S6 , SSS /OOOOOOOOOO 100, 00000000 10000, 
*0000001000000,0000100000000,0010000000000,0377777777777/ 
DATA MASK 1 ,MASK2/ 000007 7777777 ,0000000000077/ 

I F I XO ( I Y ) = IY 

MASK (IX) = IFIX0(AND(MASK1,IX) ) 

EX ( I X , N ) = IFIXO( ANDIMASK2, IARS( (N-l )*6, IABS( IX) ) ) ) 
SYMARG X , EQ ,NUM , OEN 
ATOMIC X ( 63 ) 

C INITIALIZE NM , ON , AND THE NO ARRAY 

NM = 0 
DN = 0 
K = 0 

OUT 1 = .FALSE. 

0UT2 = .FALSE. 

DO 5 1=1, M 
L = 0 

DO 1 J =1 »M1 

LET LOCATE = FIND EQ ( J ) , APP , ALL , ( X (I) ) 

I F ( .NOT. LOCATE) GO TO 1 
K = K+l 

IF(K.GT.IOO) GO TO 301 
L = L + l 

NQ ( K ) = S4+S 5* I +S6* J 

1 CONTINUE 

S = L+S2*L 
N1 = K-L+l 
N2 = K 

DO 2 J=N1,N2 

2 NO ( J ) = NQ( J )+S 
IF(NM.NE.O) GO TO 4 

LET LOCATE = FIND NUM , APP , ALL , < X (I) ) 

I F ( .NOT. LOCATE) GO TO 4 
NM = I 

4 IF(DN.NE.O) GO TO 5 

LET LOCATE = FIND DEN , A PP , ALL , ( X ( I ) ) 

IF( .NOT. LOCATE) GO TO 5 
DN = I 

5 CONTINUE 

IF(NM.EQ.O) GO TO 302 
IF(DN.EQ.O) GO TO 303 
MNO = K 

CALL SORT ( NQ, MNO ) 

NO ( MNQ+1 ) = SSS 
K = 1 

DO 9 1=1 , M 1 
L = 1 

LOCATE = .FALSE. 

6 I F ( EX ( NQ ( K ) , 6 ) *GT . I ) GO TO 7 

I F ( EX(NQ(K ) ,5) .EQ.DN) LOCATE = .TRUE. 

K = K+l 
L = L + l 
GO TO 6 

7 N1 = K-L+l 
N2 = K— 1 
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S = S3* ( L-l ) 

DO 8 J=N1,N2 
NQ( J ) = NO ( J > + S 

8 I F ( LOCATE ) NQ(J) = -NQ(J) 

9 CONTINUE 
MNEQ = 0 
MNQI = 0 

C SEARCH FOR A VARIABLE TO ELIMINATE 

100 MIN = SSS 

IF(MNEQ.GT. 146 ) GO TO 304 
DO 10 1=1, MNO 
S = EX ( NQ( I ) ,5) 

IF(S.E0.NM.0R.S.E0.DN.0R.S.E0.63) GO TO 10 
MSK = MASK ( NO ( I ) ) 

I F ( NO ( I ) .LT.O) MSK = MSK+S5 
IF(MIN.LE.MSK) GO TO 10 
IMIN = I 
MIN = MSK 

10 CONTINUE 

C ARE WE DONE 

IF(MIN.EQ.SSS) GO TO 25 
C ELIMINATE VARIABLE FROM EQUATIONS 

NO = E X ( NO ( IMIN) ,3) 

EX6 = EX ( NO ( IMIN) ,6) 

IFINO.LT. 2) WRITE (6,405) EX6 
EX5 = EX ( NO ( IMIN) ,5) 

MNEO = MNEQ+1 
NEQ(MNEO) = EX5 
MNEO = MNEQ+1 
NEQ(MNEO) = EX6 
MNQ2 = MNO 
DO 11 1=1, MNO 

11 I F ( EX ( NO ( I ),6).EQ.EX6) GO TO 12 

12 N 1 = 1-1 

NQ< IMIN) = SSS 
SS = EX6*S6 
DO 14 1=1, MNO 

IF(EX(NO(I),5).NE.EX5) GO TO 14 
IEO = EX ( NO ( I ) , 6 ) 

MNEO = MNEQ+1 
NEQ(MNEO) = IEO 
LOCATE = NO (I). LT.O 
NO ( I ) = SSS 
MNO 1 = MNQ1+1 
S = S6* I EQ-SS 
DO 13 J=1 ,N0 
L = Nl+J 

IF(NO(L).GE.SSS) GO TO 13 

MNQ2 = MN02+1 

NO ( MN02 ) = IABS(NQ(L) )+S 

I F ( LOCATE ) NQIMN02) = -NO ( MNQ2 ) 

13 CONTINUE 

14 CONTINUE 
MNEO = MNEO+1 
NEO(MNEO) = 0 
DO 15 J= 1 , NO 

L = Nl+J 
MNQI = MNQ1+1 

15 NO ( L ) = SSS 

CALL SORT ( NQ,MNQ2 ) 
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MNQ = MNQ2-MNQ1 
MN01 = 0 
K = 1 

C UPDATE S3,S4 AND COMBINE LIKE VARIABLES 

DO 21 1=1, Ml 
L = 1 

16 I F ( EX ( NO (K),6),GT«I) GO TO 17 
K = K+l 

L = L + l 
GO TO 16 

17 IF(L.EQ.l) GO TO 21 
N1 = K-L+l 

N2 = K— 1 
N = 0 

DO 19 J=N1,N2 

I F ( NO ( J ) .GE.SSS) GO TO 19 
EX5 = EX( NO ( J ) , 5 ) 

EX4 = EX ( NO ( J ) ,4 ) 

DO 18 L=N1,N2 
IFIL.EQ.J ) GO TO 18 
IF(EX(NQ(L) ,5>.NE.EX5) GO TO 18 
S = EX ( NQ { L ) , 4 ) 

IF(S+EX4.LE.63> GO TO 171 
S = 63-EX4 

IF( .NOT.OUT1 ) WRITE (6,406) EX5,I 
01JT1 = .TRUE. 

171 S = S*S4 

I F ( NO ( J ) . LT .0 ) S = -S 

NO ( J ) = NO ( J ) +S 

NOIL) = SSS 

N = N+l 

MNQ 1 = MNQ1+1 

18 CONTINUE 

19 CONTINUE 

S = <N2-N1+1-N)*S3 
DO 20 J=N1 ,N2 
SS = E X ( NQ ( J ) , 3 ) 

IFISS.E0.63) GO TO 20 
SS = S— SS*S3 

IF(NQIJ).LT.O) SS = -SS 
NO ( J ) = NO ( J )+SS 

20 CONTINUE 

21 CONTINUE 

C UPDATE SI AND S2 

DO 24 1=1, M 
J = 0 
K = 0 

DO 22 L= 1 , MNO 

I F ( EX ( NO (L),5).NE.I) GO TO 22 
J = J + l 

K = K+EX ( NQ ( L ) ,4 ) 

22 CONTINUE 
IF(J.EQ.O) GO TO 24 
IFIK.LE.63) GO TO 221 
K — 6 3 

I F ( .N0T.0UT2 ) WRITE (6,407) I 
0UT2 = .TRUE. 

221 S = J*S2+K 

DO 23 L=1 ,MNO 

IF(EX(N0(L),5).NE.I) GO TO 23 
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SS = S-S2*EX(NQ<L> ,2)-EX(NQ(L> ,1) 

IF(NQ(L) .LT.O) SS = -SS 
NQ( L ) = NO( L )+SS 

23 CONTINUE 

24 CONTINUE 

CALL SORTING), MNQ) 

MNQ = MNQ-MNQ1 
MNQ1 = 0 
GO TO 100 
C RETURN 

25 MNEO = MNEQ+1 
NEO(MNEO) = 0 
DO 26 1=1, MNO 

26 I F ( NO (I).LT.SSS) GO TO 27 
GO TO 308 

27 EX6 = EX(NQ( I) ,6) 

MNEO = MNEO+1 
NEQ(MNEO) = EX6 
DO 28 1=1, MNO 

I F ( NO ( I ) .GE.SSS) GO TO 28 
I F < EX ( NO ( I) ,6) .NE.EX6) GO TO 29 

28 CONTINUE 
RETURN 

29 WRITE (6,409) 

RETURN 

C ERROR 

301 WRITE (6,401) 

GO TO 500 

302 WRITE (6,402) 

GO TO 500 

303 WRITE (6,403) 

GO TO 500 

304 WRITE (6,404) 

GO TO 500 

308 WRITE (6,408) 

500 STOP 

401 FORMAT ( 48H **ERROR** THE TOTAL NUMBER OF TERMS EXCEEDS 100) 

402 FORMAT (49H TERROR** THE OUTPUT VARIABLE COULD NOT BE FOUND) 

403 FORMAT (48H *$ERROR** THE INPUT VARIABLE COULD NOT BE FOUND) 

404 FORMAT (27H **ERROR** SYSTEM TOO LARGE) 

405 FORMAT ( 39H Y*WARNING** ONLY ONE TERM IN EQUATION ,12) 

406 FORMAT (45H **WARNING** THE NUMBER OF TERMS INVOLVING X(,I2,14H) I 
*N EQUATION ,I2,12H EXCEEDED 63,/,12X,36H SOLUTION MAY NOT BE IN MI 
*NI MUM FORM) 

407 FORMAT (45H **WARNING*# THE NUMBER OF TERMS INVOLVING X(,I2,30H) I 
*N ALL EQUATIONS EXCEEDED 63,/,12X,36H SOLUTION MAY NOT BE IN MINIM 
*UM FORM) 

408 FORMAT ( 34H **ERROR** SYSTEM CANNOT BE SOLVED) 

409 FORMAT (45H **WARNING#* SOLUTION BELOW MAY NOT BE UNIQUE) 

END 
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SORT 


SI8FTC SORT. 

SUBROUTINE SORT(IX,N) 

DIMENSION I X ( 1 » 

M = N-I 
DO 1 1=1, M 
K = 1 + 1 
DO 1 J=K,N 
IA = IX(I» 

I F ( IABSt IAI.LE.IABSI I X ( J) ) ) GO TO 1 

mil = ix(j) 

I X ( J ) = IA 
1 CONTINUE 
RETURN 
END 


COS 


$ I BFTC TAKOUT 

SUBROUTINE COS(X) 
ENTRY S IN ( X ) 

ENTRY TANH(X) 
ENTRY AT AN ( X ) 
ENTRY ALOG(X) 
ENTRY DUMP ( X ) 
ENTRY PDUMP(X) 
ENTRY FMCDIF(X) 
ENTRY EXPDMP(X) 
ENTRY FMCDMP(X) 
STOP 
END 


NBR 


$ I BFTC AIDI 

LOGICAL FUNCTION NBR ( X , L ) 

REAL X(l) 

N = LGR(30,FC(X,L) > 

NBR=. FALSE. 

I F ( (N.GE.O) .AND.tN.LE.9) ) NBR=.TRUE. 

RETURN 

END 
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REMOVE 


SIBFTC A I D2 

SUBROUTINE REMOVE(X) 

INTEGER X < 1 ) 

LOGICAL NBR 

DATA AAST,AMINUS, AE, ALP, ARP, ADP, AZ ER0/6H*00000 , 6H-00000, 
1 6H£00000,6H< OOOOO, 6H) 00000 ,6H. 00000 , 6H000000/ ,MSW/ 1 / 
DATA ADLR,K/6H$00000,1/ 

LMAX = X(X)+6 
L = 7 

109 GO TO (110, 111, 112, 113), K 

110 IF(FC(X,L).NE.AAST) GO TO 201 
GO TO 202 

111 IF(FC(X,L).NE.AAST) GO TO 201 
GO TO 203 

112 IF(FC(X,L).EO.ALP) GO TO 204 
120 I F ( NBR ( X , L ) ) GO TO 140 

GO TO 201 

113 IF(FC< X,L> .NE.AMINUS) GO TO 201 
L = L+ 1 

MSW=2 
GO TO 120 
140 L=L+1 

I F ( NBR ( X , L ) ) GO TO 140 

114 IF(FC(X,L>.NE.ADP) GO TO 201 
L = L+ 1 

115 IF(FC(X,L).NE.AZERO) GO TO 201 
L = L+ 1 

116 GO TO ( 150,1 60 I , MSW 

150 IFU.EQ.LMAX+l ) GO TO 170 

IF(NBR(X,L).OR.FC(X,L).EO.AE) GO TO 201 
GO TO 170 
160 MSW=1 

IF(L.E0.LMAX + 1 ) GO TO 170 
I F( FC ( X, L ) .NE. ARP ) GO TO 201 
170 CALL DELETE! X,L-1 ) 

CALL DELETE ( X ,L-2 ) 

201 K = 1 

GO TO 209 

202 K = 2 

GO TO 209 

203 K = 3 

GO TO 209 

204 K = 4 

209 L = L+l 

IF(L.LE.LMAX) GO TO 109 

210 IF(FC(X,LMAX).NE.ADLR) RETURN 
K = 1 

CALL DELETE(X,LMAX) 

RETURN 

END 


t\j on 


Function FC 


$ I8FTC A I D3 

FUNCTION FC ( X ,L ) 

DIMENSION X(l) 

DATA MASK/0770000000000/ 

LM1=L-1 

IW=LMl/6+l 

FC=AND( ALS!6*M0D(LM1,6>,X( IW) ),MASK> 

RETURN 

END 


DELETE 


SIBFTC AIDA 

SUBROUTINE DELETE(X,L) 

DIMENSION X<L ) ,ZER0(6) ,BLANK<6) 

DATA ZERO /0007777777777, 0770077777777 , 0777700777777, 

1 0777777007777,0777777770077,0777777777700/, 

BLANK / 0600000000000, 0006000000000, 0000060000000, 
0000000600000,0000000006000,0000000000060/ 

LML=L-1 
I W=LM1 / 6+1 
I P=M0D ( LM1 , 6 ) +1 

X( IW)=0R(AND(X( IW) ,ZERO( IP) ) , BLANK! IP) ) 

RETURN 
END 











EVAL 


SIBFTC EVAL 

REAL ITMP,K 

DATA TWOPI/6. 2831853/ 

DATA NEXP, NS I N , NCOS ,NS I NH , NCOSH , NSQR T , NK , NT , NZ ,NWN/10*0/ 

COMMON /C EVAL/T(10),Z(10),WN (10), K(10), NEXP, NS IN, NCOS, NSINH, NCOSH, 

* NSORT 
COMMON /OMEGA/W, A,B,C,D 

COMMON /NONLIN/ EXPA < 5 ) , S I NA ( 5 ) ,COSA ( 5 ) , 

* EXPB(5),SINB(5),COSB(5), 

* EXPC(5),SINC(5),COSC(5), 

* EXPD( 5 ) ,SIND<5 ) ,COSD< 5) , 

* SINHA(5)»COSHA(5)»SQRTA(5)> 

* SINHB(5),COSHB(5),SQRTB(5), 

* SINHCI5 ) ,C0SHC < 5) ,SQRTC ( 5 ) , 

* SINHD(5) ,COSHD(5) 

NAMELIST /INPUT/ NEXP , NS I N , NCOS , 

EXPA, SINA,COSA, 

EXPB, SINB,COSB, 

EXPC, SINC,COSC, 

EXPO, SIND, COSD, 

NSINH, NCOSH, NSORT, 

SINHA,COSHA,SQRTA, 

SINHB,COSHB, SORTS, 

SINHC,COSHC,SQRTC, 

SINHD,COSHD, 

NK,K,NT,T,NZ,Z,NWN,WN,FSTART,FEND,DELTAF 


* 

* 

* 

* 

* 

$ 

* 

* 


1 READ (5, INPUT) 
WRITE (6,201) 
IF(NEXP.NE.O) 
*1=1 ,NEXP ) 
IF(NSIN.NE.O) 
*1=1 ,NSIN) 

I F ( NCOS . NE . 0 ) 
*1=1, NCOS) 
IF(NSINH.NE.O) 


WRITE (6,202) (I,EXPA(I),EXPB(I),EXPC(I),EXPD(I), 
WRITE (6,203) ( I ,SINA( I ) ,SINB ( I ) ,SINC( I ) ,SIND( I ) , 
WRITE ( 6 , 204 ) ( I , COSA ( I ) ,COSB( I ) ,COSC( I ) ,COSD( I ) , 


IF(NK.NE.O) 

IF(NT.NE.O) 

IF(NZ.NE.O) 

IF(NWN.NE.O) 


WRITE (6,205) ( I , S INHA ( I ) , S I NHB ( I ) , S I NHC ( I ) , 
*S I NHD ( I ) , 1 = 1 , NSINH ) 

IF(NCOSH.NE.O) WRITE (6,206) ( I , COSH A ( I ) , COSHB ( I ) , COSHC ( I ) , 

*COSH D ( I ) , 1 = 1 , NCOSH ) 

IF ( NSORT. NE.O ) WRITE (6,207) ( I , SOR T A ( I ) , SOR T B ( I ) , SOR TC ( I ) , 

*1=1, NSORT) 

WRITE (6,208) ( I , K ( I ) , I = 1 , NK ) 

WRITE (6,209) ( I , T ( I ) , I = 1 , NT ) 

WRITE (6,210) ( I ,Z( I ) , 1 = 1, NZ ) 

WRITE (6,211) ( I rWN( I ) , I = 1,NWN) 

WRITE (6,212) FSTART,FEND,DELTAF 
IF(DELTAF.EO.O.O) GO TO 1 
F = FSTART 
I = 0 

W = TWOP I *F 
CALL COEFF 
DIV = A 

IF(DIV.EQ.O.O) DIV = B 
IF(DIV.EQ.O.O) DIV = 1.0 
A = A/DIV 
B = B/DIV 
C = C/DI V 
D = D/DIV 
TEMP = C**2+D**2 
RTMP = ( A*C+B*D)/TEMP 
I TMP = (B*C-A*D)/TEMP 
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ATMP = SQRT( ( A**2+B**2 )/TEMP ) 

PTMP = ATAN2( I TMP ,RTMP ) *57. 2957795 
I = 1 + 1 

WRITE (6,213) I,FfA,B»C*D,RTMP,ITMP»ATMP,PTMP 
c - ECT AST4.pi HAT I I > *DFL T AF 


201 

202 


1C ( U i- 1 I * f • 

F = FSTART+FLOAT ( I ) *DELTAF 
IF(F.GT.FEND) GO TO 1 
GO TO 2 

FORMAT ( 1H1 »// »36H FORTRAN 


FORMAC OUTPUT) 


FORMAT 
♦5HEXPD 

203 FORMAT 
*5HSIND 

204 FORMAT 
*5HCOSD 

205 FORMAT 


10X»5HEXPB 

, 10X, 5HEXPC 

,10X, 

10X , 5HS I NB 

, 10X , 5HS I NC 

,10X, 

10X , 5HC0SB 

, 10X , 5HC0SC 

,10X, 


,/,< I5,4G15.5) ) 

, ( 1HO, 3X , 1H I , 5X , 5HSI NHA, 10X,5HSI NHB, 10X * 5HS INHC,10X, 

*5HSINHD,/, ( I5,4G15.5) ) ,,,- 

FORMAT ( 1H0 ,3X,1HI,5X, 5HC0SHA, 10X,5HC0SHB, 10X, 5HC0SHC, 10X, 


206 . . , 

*5HC0SHD,/,< I5,4G15.5) ) , . T _ K 

207 FORMAT ( 1H0, 3X , 1H I , 5X , 5HSQRT A, 10X , 5HSQRTB, 10X, 5HS0RTC , / , ( 15 , 3G15.5 

*) ) 

FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 

*G1 5 #5 f / t fJ.JAfJ.rir fUAi 

*4HIMAG, 10X,3HABS, 1 IX, 5HANGLE , / ) 

213 FORMAT (I5»9G14.5) 

END 


208 

209 

210 
211 
212 



7HDELTAF=, 
4HREAL , 10X , 



EXPF 


$ I BFTC FCTNS 

COMPLEX FUNCTION EXPFX(I) 


COMPLEX X » Y 


* 

* 

* 

* 

* 

* 

* 


COMMON /OMEGA/W,A,B,C,D 

COMMON /NONLIN/ EXP A < 5 ) , S I NA ( 5 ) , COSA < 5 ) , 
EXPB(5)»SINB(5)tCOSB(5)t 
EXPC ( 5 ) ,SINC(5),COSC(5), 
EXPO(5),SIND(5),COSD(5), 

S INHA (5 ) ,COSHA( 5) ,SQRTA<5) 
SINHBI 5) tC0SHB(5) » SQRTB ( 5) 
S I NHC (5)»C0SHC(5)»S0RTC(5) 
S IN HO (5)»C0SHD(5) 

X = CMPLX(EXPCm-EXPAm*W**2,EXPB(I)*W) 

Y = EXPO( I )*CSQRT<X) 

EXPFX = CEXP(Y) 

GO TO 1 

ENTRY SINFXI I ) 


♦ 


X = CMPLXISINC ( I )-SINA<I)*W**2-SINB(I>*W) 

Y = S IND ( I )*CSQRT(X) 

EXPFX = CSIN(Y) 

GO TO 1 

ENTRY COSFX ( I ) 

X = CMPLXICOSCI I )-COSA( I )*W#*2tCOSB( I )*W) 

Y = COSD( I ) *CSQRT ( X ) 

EXPFX = CCOS(Y) 

GO TO 1 

ENTRY SINHFX(I) 

X = CMPLXISINHCm-SINHAI I ) *W**2 , S I NHB < I >*W) 

Y = S INHO ( I ) #CSQRT ( X ) 

XI = REAL ( Y ) 

Y 1 = A I MAG ( Y ) 

EXPFX = CMPLX(SINH(X1 )*COS< YI ) ,COSH(Xl )*SIN( Y1 ) ) 
GO TO l 

ENTRY COSHFX ( I ) 

X = CMPLX(C0SHC<I)-C0SHAm*W**2,C0SHB( I )*W) 

Y = C OS HD ( I )*CSQRT(X) 

XI = REAL(Y) 

Yl = A I MAG ( Y ) 

EXPFX = CMPLX(COSH( XI )*COS 1 Yl ) ,SINH( XI )*SIN( Yl ) ) 
GO TO 1 

ENTRY SQRTFX ( I ) 

X = CMPLX(SQRTCm-SQRTA{I)*W#*2,SORTB< I )*W) 
EXPFX = CSORT(X) 

1 RETURN 
END 
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non o o 


COEFF 


$IBFTC COEFFR 

SUBROUTINE COEFF 

THIS SUBROUTINE ACCEPTS THE FORMAC OUTPUT IN THE COMPLEX RATIONAL 
FORM 

REAL I EXP, IS IN, ICOS , I S INH, I COSH, I SORT , K 
COMPLEX X,EXPFX,SINFX,COSFX,SINHFX,COSHFX,SQRTFX 
DIMENSION REXP(5)*IEXP(5)*RSQRT(5) , I SORT ( 5) , 

* RSIN(5),ISIN(5),RCOS(5),ICOS<5), 

* RSINH<5) ,ISINH(5),RCOSH<5), IC0SHI5) , 

* A(20) ,B(20) ,C (20) ,D( 20) 

COMMON /CEVAL/Tl 10) ,Z ( 10 ) , WN ( 10 ) , K < 1 0 ) , NEXP, NS I N, NCOS, NS I NH, NCOSH, 

* NSQRT 
COMMON /OMEGA/W,AA,BB,CC,DD 
IF(NEXP.EO.O) GO TO 2 
DO 1 1=1, NEXP 
X = EXPFX ( I ) 

REXPf I) = REAL ( X ) 

1 I EXP ( I ) = AIMAG(X) 

2 IF(NSIN.EQ.O) GO TO 4 
DO 3 I = 1 , NS IN 
X = S I NFX ( I ) 

RSIN(I) = REAL(X) 

3 I S I N ( I ) = A I MAG ( X ) 

4 I F ( NCOS . EO .0 ) GO TO 6 
DO 5 1=1, NCOS 
X = COSFX ( I ) 

RCOS< I ) = REAL(X) 

5 ICOS ( I ) = A I MAG ( X ) 

6 IF(NSINH.EQ.O) GO TO 8 
DO 7 1=1, NS I NH 
X = SINHFX(I) 

RSI NH ( I ) = REAL(X) 

7 I S I NH ( I ) = A I MAG ( X ) 

8 IF(NCOSH.EO.O) GO TO 10 
DO 9 1=1, NCOSH 
X = COSHFX ( I ) 

RC OSH ( I ) = REAL ( X ) 

9 ICOSHt I ) = AIMAG(X) 

10 IF(NSQRT.EO.O) GO TO 12 
DO 11 1=1, NSORT 
X = SORTFX ( I ) 

RSORTI I) = REAL(X) 

11 ISQRTf I ) = A I MAG ( X ) 

12 CONTINUE 


THE COMPLEX RATIONAL FORM OF THE FORMAC OUTPUT IS TO BE INSERTED 
AT THIS POINT 

C************************ *********************************************** 



AA 

= 0.0 


DO 

13 J=1 ,NA 

13 

AA 

= AA+ A ( J ) 


BB 

= 0.0 


DO 

14 J=1,NB 

14 

BB 

= BB+BIJ) 


CC 

= 0.0 


DO 

15 J=1,NC 

15 

CC 

= CC+C(J) 


DD 

= 0.0 


DO 

16 J=1 ,ND 

16 

DD 

= DD+DtJ) 


RETURN 


END 



COEFF 


$ I BFTC COEFFS 

SUBROUTINE COEFF 

C THIS SUBROUTINE ACCEPTS THE FORMAC OUTPUT IN THE S FORM 

REAL K 

COMPLEX N2,D2,EXPFX,SINFX,COSFX,SINHFX,COSHFX,SQRTFX,S 
COMPLEX EXPFI5) ,SINF<5) ,COSF( 5) , S I NHF ( 5 ) ,COSHF ( 5 ) ,SQRTF(5) 

COMMON /CEVAL/TUO) ,Z< 10),WN< 10) ,K( 10) , NEXP, NS I N, NCOS , NS INH, NCOSH, 
* NSQRT 

COMMON /OMEGA/W, A,B,C ,D 
S = CMPLX(0.0,W) 

IF(NEXP.EO.O) GO TO 2 
DO 1 1=1 »NEXP 

1 EXPF ( I ) = EXPFX ( I ) 

2 IF(NSIN.EO.O) GO TO 4 
DO 3 I=1,NSIN 

3 SINFI I ) = SINFXI I ) 


4 I F( NCOS . EQ.O ) GO TO 6 
DO 5 1=1, NCOS 

5 COSFI I ) = COSFX ( I ) 


6 

I F 1 

INSINH. 

» EO .0 ) 

GO 

TO 

8 


DO 

7 1 = 1, 

,NSINH 




7 

SINHFI I ) 

= SINHFXI I 

) 


8 

I F < 

NCOSH , 

, EO, 0 ) 

GO 

TO 

10 


DO 

9 1=1, 

, NCOSH 




9 

COSHF ( I ) 

= COSHFX ( I 

) 


10 

I F 1 

:nsort, 

.EQ.O) 

GO 

TO 

12 


DO 

li i = i 

L , NSQRT 




1 1 

SORTF ( I ) 

= SORTFX ( I 

) 



12 CONTINUE 

C THE S FORM OF THE FORMAC OUTPUT IS TO BE INSERTED AT THIS POINT 

C £ 4 4 # & ***##**#:*##* £ "S' * $ * A * * * s)c * >j< $ * if $$$$$$$$ ,}: 5jt # $ 

A = REAL ( N2 ) 

B = AIMAGIN2) 

C = REALID2) 

D = A I MAG ( D2 ) 

RETURN 

END 
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APPENDIX D 


METHOD OF DETERMINING ORDER 

The method of determining the order of equation reduction simulates the FORMAC 
reduction. Initially, certain data are obtained for each equation and variable. These 
data are updated at each step of the simulated reduction and are used to determine the 
variable to be eliminated. This method was chosen because the FORMAC expressions 
need be examined only once to initialize the simulated reduction. Also, the data retained 
at each step may be easily updated for the next step. 

Four pieces of information are retained at each step of the simulated reduction for 
each variable in each equation. Listed in order of importance, they are 

(A) The number of times the variable has been combined in the equation 

(B) The number of variables in the equation 

(C) The number of equations in which the variable appears 

(D) The number of times the variable has been combined in all equations 

These data are initially extracted from the FORMAC block-diagram expressions. 

The data are updated through each step of the simulated reduction by eliminating varia- 
bles and equations and by combining like variables in equations. 

The variable to eliminate and the equation used to eliminate the variable is chosen by 

3 2 

finding the smallest from the list of numbers made by forming the number (8 A + 8 B 
+ 8C + D) (where A, B, C, and D are octal numbers). This, in effect, finds the set of 
smallest A, then reduces this set by finding the intersection between this set and the set 
of smallest B, and so on for C and D. The variable and equation number associated 
with the ABCD found are then the ones desired. The variables of the transfer function 
are, of course, never eliminated. Equations that contain input variables are only used 
to eliminate a variable if there is no other choice. 

A list is formed consisting of the variable eliminated, the equation used to eliminate 
the variable, and all other equations that contain the variable. This occurs at each step 
of the simulated reduction. The list is then used for the actual FORMAC reduction. 

The data, A, B, C, and D are packed into one computer word along with the equa- 
tion and variable number. This word is tagged with a negative sign if the equation asso- 
ciated with it contains an input variable. This packing saved considerable storage space; 
however, limitations to the magnitude of A, B, C, and D, the number of equations, and 
the number of variables were introduced. These magnitudes, 1 to 30 for the number of 
equations and 1 to 63 for the other data, were considered large enough for most problems. 

To illustrate this method, the set of block-diagram equations given in table I is car- 
ried through the simulated reduction in table in. 
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TABLE m. - STEPS USED IN SIMULATED REDUCTION TO OBTAIN ORDER OF REDUCTION 


[Solution is for transfer function X^/Xg] 



Step 0 





Step 1 





Step 2 





Step 3 





Step 4 




Equa- 

Vari- 

A 

B 

C 

D 

Equa- 

Vari- A 

B 

C 

D 

Equa- 

Vari- 

A 

B 

C 

D 

Equa- 

Vari- A 

B 

C 

D 

Equa- 

Vari- A 

B 

C 

D 

tion 

able 





tion 

able 




tion 

able 





tion 

able 




tion 

able 




1 

1 

1 

3 

2 

2 

1 

1 1 

3 

2 

2 

1 

1 

1 

3 

2 

2 

1 

1 1 

3 

2 

2 

2 

1 2 

2 

1 

2 

1 

2 

1 

3 

2 

2 

1 

2 1 

3 

2 

2 

1 

2 

1 

3 

2 

2 

1 

1 2 1 

3 

2 

2 

2 

2 2 

2 

1 

2 


1 6 1322 1 6 1322 1 6 1322 1 |6 1322 

2 2 1222 2 2 1 2 22 2 2 1222 2*11322 

2 3 12222 3 12222 3 12222 2 1322 

i : 

3 4 1222 4 3 1322 4 1 2 3 4 5 1 1 ' 3 i 2 2 2 |6 1322 

3 6 1222 4 5 1322 

4 3 1322 4 6 1322 

4 4 1322 5 1 1222 

4 5 1322 5 5 1222 

5 1 12 2 2 

5 5 1 2 2:2 ' 

Variable to eliminate 4 5 3 I g 




Equation to eliminate 
variable 


3 


5 


4 


1 
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